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The Navy/NASA Engine Program (NNEP89) - A User s Manual 


By Robert M. Plencner and Christopher A. Snyder 

National Aeronautics and Space Administration 
I-ewis Research Center 
Cleveland, Ohio 44136 


SUMMARY 

An engine simulation computer code called NNEP89 has been written to perform one 
dimensional steady state thermodynamic analysis of turbine engine cycles. By using a very 
flexible method of input, a set of standard components are connected at execution time to 
simulate almost any turbine engine configuration that the user could contemplate. The code 
has been used to simulate a wide range of engine cycles from turboshafts and turboprops co air- 
turbo-rockets and supersonic cruise variable cycle engines . Off-design performance is 
calculated through the use of component performance maps. A chemical equilibrium model is 
rr; Q ated ^ ® det ? uate,y chemical dissociation as well as model virtually any fuel. 

V s T Sta ^ ard FO , RTRAN77 with clear structured programming and extensive 
mternal documental, on. The standard FORTRAN?? programming allows it to be installed onto 
most mainframe computers and workstations without modification, 

mvpdbo* 16 N ^ EP89 has 1)660 derived from the Navy/NASA Engine Program (NNEP). 

EF89 provides many improvements and enhancements to the original NNEP code and 
incorporates features which make it easier to use for the novice user. NNEP89 has been written 
t» execute the eld NNEP input files without ehenges to the program or the input files 

This report serves as a comprehensive user's guide for the NNEP89 code It 

TI*? 8 g ? eral de8CT j pti0 ° ofthe 6 <>de. comprehensive input description, program flow 
charts and sample input and output cases. 
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1*0 INTRODUCTION 

DeveloDm Res ** rch Cente . r in function with the Naval Air 

Engine Program, NNEP^reM)*™ code ? alled the ^^ASA 

existing Navy code NEPCOMP MnvT '!?? expanded upon the capabilities of an 

engines fo?^ d n^h^ a Zfu* material j c apabilities and recent studies of air breathing 

that NNEP could not toddle adequ^telyTim u f" 8 * 11 * CVC,es and engine editions 

these engine cvcles At y V ®P h,gh tem P*ratures are reached in many of 

streams can occur, which NNEP can gas 
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cooling bleS How r^ir T^Z 7-% ?r^, P *! Calc Y ,at f d »« nn-ount of turbine 
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Th>» document provides a general description of the NNEP89 code anrl »ll th* 
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2.0 GENERAL CODE DESCRIPTION 


The NNEP89 engine simulation computer code performs one dimensional, steady state, 
thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a 
set of standard components are connected at execution time to simulate almost any turbine 
engine configuration that the user could contemplate. Off-design performance is calculated 
through the use of component performance maps. The compressor and turbine performance 
maps are scaled by the code to match the design point pressure ratio, corrected weight flow and 
efficiency of the engine being modeled. Engines that change configurations over various 
portions of their flight regimes are modeled using a feature of the code that allows the user to 
defin 3 multiple configurations of the same engine. These engines with multiple configurations 
are called multimode engines. The default thermodynamic routine used in the code is preset for 
mixtures of air and JP4 fuel. A chemical equilibrium model is incorporated as an option to 
adequately predict thermodynamic properties when chemical dissociation occurs as well as 
when using virtually any fuel. 

In general, two separate input files are required to run the code. The first contains 
inputs which tell the code what components will be used and how those components are 
configured to form a specific engine model. Detailed inputs for each of these components 
describe the desired component model. Also included in this input file are global inputs which 
control program input/output, execution, optimization, turbine cooling, thermodynamic property 
calculations, and installation effects calculations. The second input file contains all 
performance map tables which are generally used to model off design performance of 
components such as compressors and turbines. If only design point calculations are being 
computed the map tables file is usually not required. 

Although NNEP89 only does a thermodynamic analysis, it has been coupled to a number 
of other codes to extend this capability. The capability to estimate engine weights was added to 
the program (ref. 10). Presently, this is accomplished by using additional program libraries to 
perform the weight analysis. The weight program individually weighs each component, using 
user input for many design parameters.. There is an option for graphical output of the flow 
path and interactive weight determination. 

Simplified installation effects calculation are built into the program to give a 
preliminary estimate of inlet drag and nozzle drag. As an alternative to these internal 
installation calculations there is an additional pi ugram library which performs more detailed 
installation effects calculations (ref. 11). This additional program library can utilize much more 
sophisticated inlet and nozzle performance maps, to more accurately estimate installation losses 
and their effect on overall engine performance. The installation package can also perform much 
more detailed inlet and nozzle weights if sufficient inlet and nozzle system information is 
known. 
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3.0 COMPONENT PERFORMANCE MAPS 

One of NNEP89's greatest strengths is its capacity to interpolate user-supplied data 
tables to estimate component performance. These data tables consist of a dependent variable 
which is a function of up to three independent variables ( dependent variable = F (x,y,z) ). 
Consider the table data as three dimensional, composed of a series of planes with each plane 
assigned a value call Z. Then, on each Z plane, the dependent variable (ordinate axis), is a two 
dimensional function of X (abscissa axis) and Y (see Figure 1). These three dimensional tables 
are often referred to as maps. The data format for inputting all tables is given in Section 10.3. 



Figure 1. Three Dimensional Map Representation 


X 


• Yi 



Most of the engine components have inputs which may be specified in terms of a user 
input tabular function of one to three independent variables. The specific independent variables 
used are predetermined by the program. These functions are discussed in the descriptions of 
each of the components in Section 6. The user can also set up a variable schedule, which allows 
the user to set up a data table which sets any input variable to be a function of up to three other 
component inputs, component outputs, or engine flow station properties that the user specifies 
through the van ’ * schedule inputs. 



All the map bles required for a given engine model must be contained in a separate 
map file. This entire map file is read in by the program at the beginning of the execution of a 
set of cases and the entire file is stored in an array. This array is later interpolated by the code 
using a cubic spline interpolation technique to get the values of the depe ndent variables 
corresponding to the current values of the independent variables. If ar. independent variable is 
outside the range of the tabular data, a linear extrapolation will be performed using the slope of 
the spline at the last point in the table. 
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pressure ratio, the mass flow through this turbine, and user specified design point on the map, 
the program calculates the turbine map scale factors. The map scale factors for the turbine 
component have the same definition as the compressor map scale factors. These map scale 
factors will be used by the program for off-design points to convert the turbine corrected mass 
flow and efficiency interpolated from the performance maps into the actual performance 
properties used by the program. 

There is a separate program available for generating turbine performance maps (ref. 
14). The program is easy to use and can generate turbine performance maps with operating 
characteristics tailored for the user’s application. These generated performance maps are 
already in the NNEP89 format. 


Many of the engine component models are already set up to use performance maps for 
certain input variables. (For example, the inlet recovery data table is set up to use a table of 
inlet total pressure recovery as a function of inlet referred flow and Mach number ) The format 
of these data tables is the same as any other NNEP89 data table input, such as the compressor 
and turbine maps. The user can also input a user-designed data table and use the variable 
schedule program control component to have the program set any input specification of any 
engine component as a function of any 3 component inputs, outputs, or engine performance 
properties. An example how to do this is include in the Section 10.2.16. 


12 


w 




INPUT FIELD 

DESCRIPTION 

USAGE 

TITLE LINE 

1 line input identification 

Always required 

&D &END 

Global &D Namelist Inputs 

Always required 

REACTANTS 

CEC switch to read in reactant 
data 

Only input when using the 
CEC option 

(CEC reactants) 

Reactant data information 

Only input when using the 
CEC option 

(blank line) 

Signals end of CEC reactant 
data 

Only input when using the 
CEC option 

THERMO 

CEC switch to read in thermo 
data 

Only input when using the 
CEC option 

(CEC THERMO DATA) 

Thermo data information 

Only input when using the 
CEC option 

END 

Signals end of CEC thermo data 

Only input when using the 
CEC option 

NAMELIST 

CEC switch to read in CEC 
Namelist data 

Only input when using the 
CEC option 

&INPT2 &END 

CEC Namelist data 

Only input when using the 
CEC option 

END (CEC switch to signify 
end of CEC input data) 

Signals end of CEC input data 

Only input when using the 
CEC option 

&D MODE=l &END 

SPEC and KONFIG Namelist 
inputs for MODE=l 

Always required 

&D MODE=2 &END 

SPEC and KONFIG Namelist 
inputs for MODE=2 

Only input for engines with 
2 or more modes 

&D MODE=3 &END 

SPEC and KONFIG Namelist 
inputs for MODE=3 

Only input for engines with 
3 or more modes 

&D MODE=4 &END 

SPEC and KONFIG Namelist 
inputs for MODE=4 

Only input for engines with 
4 or more modes 

&D MODE=5 &END 

SPEC and KONFIG Namelist 
inputs for MODELS 

Only input for engines with 
5 or more modes 

&D MODE=6 &END 

SPEC and KONFIG Namelist 
inputs for MODE=6 

Only input for engines with 
6 modes 

&D &END 

1 

First off-design case or new 
design point case (these may be 
repeated as many times as 
required 

Required for off-design 
cases or multiple design 
point cases (set NCODE=3 
for multiple design points) 


Table 1. Example Input Data Format 
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5.0 NNEP89 ENGINE MODEL DEVELOPMENT * 

t 

5.1 CONFIGURING AN ENGINE 5 

Tlie NNEP89 computer code does not contain any preset engine cycles, such as a single 
spool turbojet, that the user can simply invoke through a single input. Instead, it requires the 

user to identify the mechanical and thermodynamic connection between engine components ‘ ' 

through a set of inputs. In addition, at off-design point operating conditions, the user must J 

identify potential sources of flow mismatch, work imbalance, and rotation^ speed mismatch 
within the engine through additional inputs. These mismatches and unbalances are commonly 
referred to as "errors ". Finally, the free variables that will be used to eliminate these "errors” 
must be specified by the user. 

5.1.1 Components and Flow Stations 

To use the code, the user must specify all the components which are contained in the 
desired engine configuration. Unique component numbers are assigned to each of these 
components. Since the components do not have to be labeled sequentially, flow stations are 
defined for flows entering and leaving each component. These flow stations tell the code how 
the components are connected together; the downstream flow station number of one component 
must correspond to the upstream flow station number of the component that immediately 
follows it. The restrictions that apply to the assignment of component numbers and flow 
stations are given in Table 2. 


The primary inlet must be always be given a component number of one. 

The primary upstream flow station of this inlet must be one. 

Component numbers must be in the range of 1 - 199. 

Flow station numbers must be in the range of 1 - 99, 

Table 2. Input Labeling Restrictions 

The thermodynamic properties of the flow are calculated at each of the flow stations 
through the engine. Table 3 lists these thermodynamic properties and the corresponding 
Station Property Number which is used to refer to a specific property at any given flow station. 
Notice that Mach number and static pressure are stored in Station Properties 6 and 7, 
respectively. However, in general, areas are not calculated by the code. Therefore, Mach 
number and static pressure are only calculated for a few components. Mach number and static 
pressure are calculated for the flow stations entering and leaving a mixer. The nozzle 
component calculates the static conditions for ihe throat and exit flow stations. The nozzle 
throat Mach number and static pressure are stored in the flow station location corresponding to 
the nozzle entrance. The duct (or burner component) has an option to calculate the static 
conditions for the flow stations entering and leaving that component. The free stream Mach 
number into the inlet is stored as Station Property 6; however, the free stream static pressure 
and temperature entering the inlet are stored in Station Properties 2 and 3 (total pressure and 
temperature). For all other flow stations the static pressure and Mach number are not calculated 
and will show up as zero on the output. 
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Station Property Number 

Definition 

i 

Weight flow, lb/s 

2 

Total pressure, psia 

3 

Total temperature, °R 

4 

Fuel-to-air ratio 

5 

Referred flow,wVr/p Note: corrected flow, v/yje/b, is 
printed on the output W\/t/p = 1.5497 w>/e/8 

6 

Mach number 

7 

Static pressure, psia 

8 

t — 

Interface corrected flow error 


Table 3. Station Property Definitions 


To aid in assigning the required inputs, a block diagram of the initial engine 
corfiguration, as shown in Figure 5, is typically sketched by the user. The component numbers 
and flow stations are added to this block diagram to aid in visualization of the engine layout 
and to help in formulation of the input file. The component numbers and flow station numbers 
in this example do not follow any logical sequence to illustrate that they can be assigned in an 
arbitrary manner (subject to the restrictions in Table 2). However, it is recommended that 
when possible the component numbers and flow station numbers be labeled sequentially to aid 
the user in reading and understanding the results. 


♦ 
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Figure 5. Sample Two-Spool Turbojet 
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configuration of the desired engine. Hie KONFIG array location definitions are defined for 
each component type in Section 10.2. The detailed design information required to define each 
component is input through the Namelist variable SPEC. SPEC is a two dimensional array of 
15 by 200. Each of the components has 15 design input parameters associated with it. 11)686 
detailed definitions are specified for each component in Section 10.2. 


Component Type 

Quoted string 
Input 

I Engine Components 

Inlet 

TNLT 

Duct or Burner 

DUCT 

Water injector 

WINJ' 

Gas Generator 

'GGEN' 

Compressor or Fan 

COMP' 

Turbine 

TURB' 

Heat Exchanger 

HTEX’ 

Flow Splitter 

SPLT 

Flow Mixer or Ejector 

M1XR' 

Nozzle 

NOZZ* 

Load or Propeller 

LOAD' 

Shaft 

SHFT 

I Program Control Components 

Variable Control 

] 'CNTL' 

Variable Optimization 

OPTV 

Variable Limit 

'LIMV 

Variable Scheduling 

’SKED' 

Conditional Control 

'VCNT 


Alternative 
Numeric Incut 



Table 4. Input Component Name Definitions 


When configuring an engine that requires a bypass or bleed flow stream, a compressor, 
duct (not a burner), or flow splitter may be used to create this secondary flow. To create this 
secondary flow stream the flow station number of the secondary stream is simply input in the 
fifth array location of the KONFIG variable for that particular component. A secondary flow 
that is bled off a compressor or a duct must either be bled overboard or it must reenter the main 
flow stream in a turbine or another duct. A secondary flow that is created by a duct or 
compressor may not pass through any other components between the point where it is bled and 
the point where it reenters the engine or is bleed overboard. If the flow is recombined with the 
main flow stream, the downstream flow station of the compressor or duct where the secondary 
flow stream is created must be the secondary upstream flow station of the duct or turbine where 
the flow streams are being recombined. However, all the flow does not have to reentev the 
engine at one location. It may reenter at multiple ducts and turbines. The momentum of the 
secondary stream is not accounted for when it is mixed in the duct or turbine with the main 
flow stream. The momentum of the secondary stream is simply set to the same momentum as 
the main flow stream at the point they are recombined. However, the enthalpy of the secondary 
flow stream is accounted for. The temperature of the recombined flow will be determined by the 
enthalpy of the two streams that are being combined. If overboard bleed is desired, the 
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secondary flow station is not connected to a duct or a turbine. This flow will not enter into any 
other calculatir ns and is simply "lost" by the code. y 

.... If ** v US€r d ** ire8 40 ha y e other components in the secondary flow stream, a flow 

r 06 USed ? T ate t S e secondary flow 8tre am. AH secondary flow streams that are 
ge erated by means of a flow splitter normally will be recombined with the main flow stream in 
a mixer or an ejector or will exit the engine in a separate nozzle. The secondary flow from a 

fl ? W 8tream totough a duct or turbine; however, the 
momentum of the secondary stream is not accounted for as was discussed above. The secondary 
flow created by a splitter must either be recombined with the main flow stream or exit toroueh* 

ItSshS’ 11 T n0t b J dUmP * • overboard Ifthe is recombined with the Lrin fW 
stream, the last downstream flow station of the components in the secondary flow stream must 

rec^b^eT d Note U th t r ai m™/ * component where th « Aow sto^ams are being 

recombined. Note that splitter/mixer combinations can be nested. In this case the furthest 

downstream splitter must reenter at the furthest upstream mixer. 

IL£ -MULTIMODE ENQTNF. CQNFT0ttt?attomq 

when . Mod *! in ? °* : multimode engines is a capability that was developed specifically for NNEP 

model T! ?° m ** T 8 *™ 1 Nav * NEp COMP code. This capability allows NNEP89 to 

model engines tha . change configuration over various portions of their flight regimes An 

example of a multimode engine is a convertible engine for high speed nSTn toe tokeoff 
rotor ° pera J e8 ® 8 a turboshaft -tli power being extracted to drive a helicopter type 

i » m °? e f engine 0perates a8 a turbofan disconnecting the helicopter 

rotor and connecting a fan stage to provide forward thrust.This would be modeled by defining 

renr^tto l ? eS S?*f d » bel °«' The fir8t m0de WOuld include a11 ^ components necessary to 
f 16 c o n flguration and a second mode would include all the components 

necessary to model the turbofan configuration. Most components would be common to both 







j . Wh ® n 8 ^ ti 1 n . g " p a multimode engine, a separate Namelist is used for each mode 
i'S: f c Tab,e 1 as well as Sample 5 (Section 12 . 5 ) for examples of h^to seHp the 
P . in ^ C 0 T ne l WhlCh 1S common 40 two modes must have the same component 
mnd^fh *^ th m0deS ' 11,8 user 18 cautioned that if a component appears in more than one 

a^imed by timoroeram 1 S* f ^ Va j aWe in * e last ™ de * ** input will be yZlTtoat 

than twn Pmgram. Up to six modes may be defined in an engine model, although more 

MODESNv^riabTe 6 3?* requ,red , The P™ 4 configuration is defined using toe 

sa by MlUni v,ri ” b, ‘ mode — - «- 

5 ^. OFF-DESIGN OPEWATTr> N 


pnniw 0 F n °. r de8i f P? int calculations NNEP89 will automatically ensure continuity of mass 
rgya d speed between components. For off-design operating points toe continuity of m’aa 
and Energ y are n.t automatically met. The uur m^determin^S “^7^" 


mass 
a 


Tim con^^ ytf Bf7*Ev»ry mm7m»nt n!!— ' 

compon«nt map and war, componant that choke, will only pa„ . ^ , monnt 
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flow, w>/§/8. For a compressor map the independent variables corrected speed, nA/5. and R 
uniquely determine the corrected flow through the compressor. However, this does not 
automatically match the flow being passed to it by the component immediately upstream from 
it. Similarly, the corrected flow of the turbine will be determined by its pressure ratio and its 
corrected speed. A nozzle with a fixed throat area will limit the corrected flow to the amount it 
can pass with Mach one at its throat. Therefore, the entrance to every compressor, turbine and 
fixed area nozzle will have the possibility of a flow mismatch. The flow coming out of the 
previous component may not be equal to the amount J flow that the component will pass. 

These imbalances are commonly called "flow errors". 

The second kind of possible mismatch within the engine is a work imbalance. The sum 
of the energy taken out of the shaft by compressing the air or by driving external loads or 
propellers may not equal the energy being put into the shaft by the turbines. This type of 
imbalance is commonly called a "work error". All shafts have potential work errors. 

Two components have special operating conditions that introduce unique errors within 
that component. For the mixer a control component is usually required to force the static 
pressure of the primary stream to equal the static pressure of the secondary stream. To design 
a heat exchanger, the input value of temperature rise must be varied until it matches the 
calculated temperature rise at the design point conditions. 

Once all the mismatches or errors are found within the engine, the user must pick an 
equal number of free variables within the engine. These free variables will be varied by the 
program in order to obtain error values that are smaller than the tolerance specified by the 
user. The free variables will depend on the engine configuration, and the particular set of free 
variables will somewhat affect the operating conditions predicted by the code. A typical list of 
possible free variables for each component type is given in Table 5. 

The variable control component is used to input the engine errors and the free variables 
into NNEP89. Determining the engine errors, choosing a set of free variables to be used and 
writing the variable control component inputs is one of the more difficult tasks when using 
NNEP89, especially for a new or infrequent user. As a result, an automatic control formulation 
option is available in NNEP89. This option is invoked by setting ACTL to be one or two. The 
program will then determine the engine errors and will choose a set of free variables to 
eliminate these errors. All the errors discussed above will be controlled except for the mixer 
static pressure balance. Because the operation of the mixer is very dependent upon the desired 
use, the user must input any desired control component for the mixer. Even though the 
program will assign all the other needed control components, the user should monitor what free 
variables are being used to ensure the engine is being modeled as the user desires. When 
executing with this option turned on the code will not try to use free variables that have been 
previously used in a user written control, nor will it try to eliminate an error that the user has 
already written a control to eliminate. Thus even if the controls are generally input by the user, 
the automatic control option may be invoked to ensure that an engine mismatch isn’t missed 
through a user oversight. 

The above discussion has shown the necessity of using the variable control component to 
eliminate mismatches within the engine. However, the use of the control component is much 
broader than that. It may be used to obtain a specified value for any input or output variable by 
varying an independent free variable within the engine. One example of this is to force the 
engine to operate at a specified compressor surge margin. Another example is to obtain a 
desired thrust level by varying the engine airflow or the burner outlet temperature. These 
controls may be used at design point, at off-design or both. Variable control components are 
used extensively in modeling engines with the NNEP89 code. 



m 


Component Tvoe 

Possible Free Variable 

Inlet 

Airflow 

Compressor 
or Fan 

R Value 


Third Dimensional Value on Stacked Maps 

Turbine 

Pressure Ratio 


Third Dimensional Value on Stacked Maps 

Burner 

Burner Outlet Temperature 

Mixer or Ejector 

Primary or Secondary Flow Area (if 
variable) 

Nozzle 

Throat Area 

Shaft 

Shaft Rotational Speed 

Heat Exchanger 

Primary Flow Temperature Change 

Splitter 

Bypass Ratio 

Load 

Horsepower Extracted 


Table 5. Possible Free Variables 


r,nd * “ ,utiM 

r P :p t s ^ 

has been reached the code will nrinf « wnen the maximum number of lteratioi 

next case. 1,1 Pnnt 8 W8rmng messa ^ e the and then proceed to the 
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i " dividu “ 1 com P on ents during ofT-design 
If IDONE(N)=0 component Nhafn^ ^ P gra . m ,f f con »ponent has been designed or not. 
for each component. For off-design points thrus^SttlroNE^^yS^ * 10 T 

z£r sss. Zt/ Z::tz'X1 on pre ” nt *^^£5^ 

switches modes, but thif new mode j” ° ne particu,ar ">«*•. user 

designed. The program will set IDONE(N)=0 for thiYnozzfe that has not yet been 

case based on current conditions ^ . le and l w,n redesi & n during that 

directly by setting IDONE(N)=l. ’ €r Can ,nput the nozzle design parameters 
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Setting IDONEsO for every component will cause the program to redesign all currentl 
used components. This has the same effect as setting IWAY=1. IWAY is discussed further in " 
Section 10.1.2. 

5-5 CONFIGURATION LIMITATIONS 

Presently, NNEP89 is capable of calculating the steady-state design and off-design 
thermodynamic performance of engine cycles with as many as 100 flow stations and as many as 
200 components. The maximum of 200 components includes engine components and program 
control components. The user is limited to using a maximum of 50 components of any one type 
(e.g. maximum of 50 ducts). Note that only 10 variable optimization components may be active 
at any one time. A maximum of 70 user-supplied component performance maps, with a 
maximum of 20,000 values may be used. 
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6.0 ENGINE COMPONENTS 


. . . . ^ Piously stated, NNEP89 has 11 engine components. The engine components are: 
inlet, duct or burner, water injector or gas generator, compressor, turbine, heat exchanger 
splitter, mixer or ejector, nozzle, load or propeller, and shaft. These components will now be 
d ti ai ; defmitlon of each in P ut variable for these components is given in 

i Ki r' M ° St 0 / the c° m Ponents contain inputs that are either a fixed variable value or a 
^ol C . member. If a table reference number is given, the program will interpolate the 

* K1 8 T, hed ta i ble US i n ? the current values of the independent variables of that particular 
table. A description of these table formats (the dependent and the independent variables) is 

f-Wr d ? scu f sions ; ^ must supply the table which relates the program 

specified dependent and independent variables in a separate map file which gets read in on unit 


&1 -INLET 

•n captures the airflow that enters the engine. The inlet conditions can be 

specified either as temperature and pressure of the airflow or as the geometric or geopotential 

^ aYmnsnh Um - Z 7°°^ U&]titude is in P ut > this component will use its internal 
m °?w determine entrance conditions. There is an option to add a 
delta to the atmospheric model temperature to simulate a non ideal day. The inlet component 
isentropically compresses the inlet flow to stagnation conditions and then applies a pressure loss 
using an inlet total pressure recovery factor. The inlet recovery may be input by the user as a 
wnstant value or as a table as a function of inlet exit corrected flow and Mach number. If no 
recovery is input, a default inlet pressure recovery schedule is used as follows: 
for Mach < 1 Recovery =1.0 

for Mach > 1 Recovery = 1 - 0.075 (Mach - l) 1 - 35 

(ib/S) C f n 7 S P ecified as either corrected airflow at inlet exit, corrected airflow at 
mlet entrance or actual airflow. The program is preset for two kinds of inlet performance maps 

inle/eritrefe^edT ‘"S **** d ™ de . d . by dynamic Pressure as a function of Mach number and 
Mach 4, second is inlet recovery as a function of inlet exit referred flow and 

the coir. Sra^ToSuon * ** C0 " ,, ’ 0ner,t "’ a ‘ > taMe referenre "'“" bCT 

6.2 DUCT o r BURNER 

A duct can be used as an engine's main combustor, an afterburner, a passage to move 

rad rSrr.d a d1*11 “ 10 an0 . thCT " * If™' *° preMurc 1 » s * 8 combustor 
and the afterburner add fuel to a gas stream to change its temperature. The fuel heating value 

chenSl bU TK effiC1 !K Cy be input * * burner is desired - unless the user is using the 1 

*^ ,bnUm the rmodynamics model. The desired combustion temperature nfay be 
J component will find the fuel to air ratio, or the fuel to air ratio may be 

of thfafr to d li he K c ® m P° nent wil1 find toe temperature. There are also terms to allow some 
. S . . be unburned, (e g. for combustor wall cooling). This unburned air will be added 
back to the main stream before the burner exit, such that the burner exit temperature will not 

known^the h'V* 1 bu ™ er te ™ peratu ™- If toe burner cross-sectional area or Mach number is 

•Sdra to STJSSJ Ea,lei8h totil «“ h “«"* * «•» 

ki j tr^n “ P ? 8Sag ®' a duct can 1)6 U8ed afl a transition from one component to the next to 
bleed off flow for turbine bleed, or bypass flow around other components Shout u .?ng a 

splitter-mixer combination. When using the duct to bypass flow, the program does not account 
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for the pressure differences from an entering bypass stream and the main stream. Total 
enthalpy will be conserved, but the mixed stream will be at the main stream’s pressure. The 
program does not follow the flow path of a duct bypass. Thus, there can not be any other 
components in-between where the bypass exits the duct and the bypass stream re-enters either 
a turbine or another duct. If the user wants the program to follow the flow path of a bypass 
stream (for example to have a bypass go through any one or several component before re- 
entering in another duct) the splitter components must be used. 


1 


The program is preset for three different kinds of user-input duct or burner performance 
maps; duct or burner pressure drop as a function of the entrance corrected airflow, combustion 
efficiency as a function of entrance referred airflow and fuel-to-air ratio, and fuel heating value 
as a function of desired temperature and entrance pressure. These maps are used by inputting 
the appropriate component map table reference number in the correct SPEC array location. 

6.3 GAS GENERATOR (Only With Equilibrium Option ON) 

The gas generator component is modeled as a rocket sending its exhaust into an engine. 
It can only be used with the chemical equilibrium option turned on. The gas generator is used 
to: 


1) Add a certain amount of a fuel or a fuel and oxidizer stream to an incoming flow, and then 
determine the total mixture's temperature, 


2) Add a fuel stream at a certain temperature into an engine, or 


® /rif s P ec ^ e< ^ an< ^ oxidizer streams together for a desired temperature or mixture ratio 
(like a rocket) into the user's engine. If the user specifies the temperature, the component 
determines the fuel-to-oxidizer mixture ratio, if mixture ratio is specified, the component will 
determine the mixture's equilibrium temperature. 


The incoming flow stream determines the gas generator’s pressure. If the upstream 
mass flow rate is zero, the user may input the pressure. The user must specify the fuel’s 
composition using the FARRAY array (discussed in Section 10.4). If there is also an oxidizer 
flow, the user must also specify the oxidizer composition using the OARRAY array (discussed in 
Section 10.4) Otherwise, the program will assume that there is no oxidizer. Gas generators are 
commonly added to a cycle by using the splitter component, with the bypass stream as the gas 
generator input stream. To model a rocket with no upstream flow the user must specify a zero 
bypass on the splitter. There are no preset performance maps for the gas generator component. 

fi.1 WATER INJECTOR (Only With Equilibrium Option OFF) 


The water injector component is used to estimate the effects of water injection into the 
engine. It can only be used with the chemical equilibrium option OFF. The specific heat, gas 
constant and the ratio of specific heats for the mixture are accurately accounted for (neglecting 
chemical dissociation) in the water injector component using a special water-air thermodynamic 
table. However, downsteam the water is accounted as a fuel. The thermodynamic properties 
are calculated based on the equivalent fuel-to-air ratio. The user can input the amount of water 
injected (up to 10%, by mass, of the gas flow entering the water injector) and the fraction of 
water that is vaporized or the user can let the program determine the amount of water 
necessary for saturation (up to 10%). The program will add the amount of water requested, and 
calculate the new gas stream temperature. There are no preset performance maps for the water 
lryector component. 
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6.5 COMPRESSOR or FAN 


The compressor component can be used for modeling both compressors and fans. At 
design point, the user specifies the desired pressure ratio, adiabatic efficiency, amount of 
compressor bleed, and how much compression the compressor bleed obtains. Compressor bleed 
can be removed at the beginning of the compression process, at the end, or somewhere in- 
between. The default assumes the bleed is removed at the exit conditions. If the user is using 
compressor performance maps, the user must tell the program the location of the design point 
on that map. The map scale factors are set at design point. Compressor performance maps and 
map scale factors are discussed in more detail in Section 3.1. For off-design points, the 
compressor component will interpolate from the compressor performance map, using the map 
scale factors to determine the compressor's operating characteristics. 

6.6 TURBINE 


At design point, the user inputs the desired adiabatic efficiency and, if there are 
multiple turbines on one shaft, the fraction of total power required by the shaft that will be 
produced by this turbine. The turbine component will determine the pressure ratio required to 
produce the necessary amount of power. The user can also set the fraction of turbine cooling 
used to total bleed flow available. The program can also estimate the amount of compressor 
bleed air required for turbine cooling at those operating conditions. This option is discussed in 
reference 9. The turbine cooling flow calculation has been further enhanced by allowing the 
user to directly input the desired vane and blade temperatures. If the user is using turbine 
performance maps to predict off-design performance, the u ar must also specify the turbine 
design point on that map. At design point, the user inputs will tell the program the desired map 
design point and the program will determine the map scale factors. Turbine performance maps 
and map scale factors are discussed in detail in Section 3.2. For off-design points, the turbine 
component will interpolate from the turbine performance map, using the map scale factors to 
determine the turbine's operating characteristics. 

6.7 HEAT EXCHANGER 


A heat exchanger transfers thermal energy from one stream to a lower temperature 
stream. The user inputs the pressure losses for both streams, the heat exchanger effectiveness, 
and an initial guess the temperature rise of the main stream. The heat exchanger effectiveness 
is used with the incoming temperatures of the two streams and their relative energies (mass 
flow*specific heat) to determine the temperature rise of the main stream. This determines the 
main stream exit temperature, which determines the amount of thermal energy transferred and 
the temperature drop of the secondary stream. At design point, NNEP89 must be configured 
with a Variable Control Component which either varies the heui exchanger effectiveness or 
varies the guess for the temperature rise of the main stream in order to force the guess value of 
temperature rise to be equal to the calculated temperature rise. This is used to "design" the 
heat exchanger. This control is not necessary for off-design points. 


The user can also cool the main stream by inputting a negative temperature rise 
(effectiveness is still a positive number and the secondary stream will be heated). The user can 
also input the heat exchanger entrance thermodynamic properties for one or both streams when 
the equilibrium option is used. This could be necessary when the upstream conditions ' r one or 
both «treams are outside the range able to be calculated by the equilibrium thermodynamic 
routines. 


The program is preset for two different kinds of heat exchanger performance maps. The 

first is a pressure loss for each incoming stream as a function of its referred mass flow. The 
second is the heat exchanger effectiveness as a function of actual main gas mass flow and the 
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ratio of secondary to main mass flows. These maps are used by putting the component map 
table reference number in the correct SPEC array location. 

6.8 SPLITTER 

Hie splitter component is used to create secondary flow paths that will be followed by 
the program. (See preceding discussions of secondary flows in Section 5.1.2 for the proper use of 
the splitter component.) For these secondary flows, the program will follow the secondary 
streams through other components, keeping track of the stream’s temperature and pressure. 

Hie user inputs to the splitter are the bypass ratio (ratio of secondary exit mass flow to primary 
exit mass flow) and the total pressure loss of the primary and secondary exit streams. There 
are no preset performance maps for the splitter component. 

6.9 MIXER or EJECTOR 

This component combines two streams with possibly different pressures, temperatures, 
and mass flows. It solves for the exit conditions by simultaneously solving the one-dimensional 
equations for the conservation of mass, momentum, and energy. The solution of these equations 
gives two results, a subsonic and a supersonic solution (so named because of the state of the exit 
stream). The supersonic solution generally results in a higher exit total pressure than the 
subsc.iic solution, but it is questionable whether the supersonic solution is physically possible. 

It is therefore left to the user to choose which solution, subsonic or supersonic, the mixer 
component will calculate. Mixing losses must be input by the user and are in the form of 
momentum coefficients for the entering streams or a momentum coefficient for the exit stream. 
There are options to control mixer inlet and exit areas, mixer design inlet static conditions, and 
injection angles for the entrance streams. The options are listed, along with some discussion of 
their use, in Section 10.2.9. There are no preset performance maps for the mixer/ejet xor 
component 


The nozzle component converts the potential energy of engine exhaust gases into kinetic 
energy in order to produce useful thrust. The nozzle component can handle both converging 
and converging-diverging (C-D) nozzles. Traditionally, the converging-diverging nozzle analysis 
in NNEP and NNEPEQ continuously varied the nozzle exit area such that the flow was fully 
expanded to ambient static conditions. This method yields the maximum thrust, but it also 
implies a variable geometry nozzle. The nozzle routine also handles fixed area ratio C-D 
nozzles. When the exit to throat area ratio is fixed, the nozzle will yield maximum thrust only 
at its design pressure ratio. If a fixed area ratio nozzle is operating at a pressure ratio less than 
® r greater than its design point pressure ratio, the nozzle is said to be overexpanded or 
underexpanded, respectively. Due to the unequal exit static and ambient pressures, a pressure 
thrust or a pressure drag term must be included in the gross thrust calculations. No separation 
criteria are applied for overexpanded nozzles; hence, the pressure drag term applied may be 
unusually large. The user also has the option to specify the nozzle exit pressure. The nozzle 
will expand the nozzle flow to this desired static pressure, regardless of whether the flow is 
overexpanded or underexpanded. Note that losses due to expansion to other than ambient 
conditions are categorized as thermodynamic losses. As always, any losses due to viscous effects 
must be accounted for by velocity coefficients and/or discharge coefficients furnished by the 
user. These loss coefficients are only applied to the momentum thrust calculations. The nozzle 
component can also apply divergence (i.e., non-axial) gross thrust losses to converging-diverging 
nozzles (ref. 7). These losses are accounted for by applying an analytically derived, geometry 
dependent divergence loss coefficient to the gross momentum thrust. 
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The fixed arearatio and divergence loss calculations are invoked only when certain 
inputs are specified. Thus old input files executed with NNEP89 will produce the same results 

DerfoiIl.JncA l0US VerS10n 1 8 the The n02zle component is preset for two kinds of nozzle 
forlf tW aps; . nozzl ® fl< 7 coeffi «ent as a function of nozzle total to ambient pressure ratio 

ifWl! EE mPUtS 8 7 Ue f ° r nOZZ1 * ewt Static pressure > nozzle flow coefficient as a function 
of nozzle total pressure to user input exit static pressure ratio) and nozzle velocity coefficient as 

***** t0 A f mbient P ressure ratio and the ratio of nozzle exit aJa to throat 

SrecISy Sor US y PUUing C ° mp0nent map table reference number in the correct 

LOAD or PROPRf .T,F!P 

ln . COmP ° nent A S used add or subtract horsepower from an engine cycle. If used as a 

oad, the horsepower subtraction must be input by the user as a negative number or be defined 

Z mSP ’ Wth horse P° wer as a Unction of actual shaft RPM. These maps are used 

y 1 putting the component map table reference number in the correct SPEC array location If 
^ component is used as a propeller, propeller operating maps can be included 

reftrencT? pr ° pel,er P erforman ce maps is discussed in greater detail in 

6.12 SHAFTS 

flnH .7?® Shaft C ,? mp °^ nt is for the mechanical connections between compressors turbines 
and loads or propellers. The user must input actual shaft rotational speed gear ratios betwepn 

the'shtdEi* 8 ^ onn ® cte<1 to th e shaft and the shaft itself, and the mechanical Efficiency betoeen 
COn " ected com P° nen t. Four components may be connected to any one shaft 

iwo ° r — -»• - ■» — -Hitt*- 



7.0 PROGRAM CONTROL COMPONENTS 


As previously stated, NNEP89 has 5 program control components. The program control 
components tell the code how to execute the cycle. These components are variable control, 
variable optimization, variable limit, variable scheduling, and conditional control. These 
components will now be discussed in detail. The definition of each input variable for these 
components is given in Section 10.2. 


7.1 VARIABL E CONTROL 

Since NNEP89 does not have any preset engine configurations, the user must define as 
input to the program any flow or work errors that the program must drive to zero. The user 
may also specify that some components operate at a specified condition (such as operating 
compressors at a constant surge margin). Controls define to the program what independent 
variables to change to get the engine model to operate in the desired manner These 
components are often the most difficult to set up for the user. The philosophy of variable 
controls was previously discussed in Section 5.3. An option is available to automatically set up 
the controls at execution time. This option is discussed in detail in reference 8 and an example 
of its use is giver, ir. Section 12.2. 


11 VARIABLE OPTIMIZATION 

T ^ ie optimization component allows the user to maximize or minimize a "cost function" 
(i.e., the variable being optimized). Through the optimization component, the user tells the 
program what variable, or variables, to change to maximize or minimize one "cost function" Up 
to 10 independent variables can be used at one time. The "cost function" can be any 
performance property or component output in the program (it can not be a station property). 
Variable optimization is discussed moie fully, with an example, in Section 10.2.14. 

11. VARIABLE LIMITER 


The variable limiter component tells the program to "watch" certain variables and t°U 
the user if the variables exceed the limits set by the user. If the user is executing with 
optimization and a limit has been exceeded, the program will penalize the optimization "cost 
function (discussed in Section 7.2) to try to get that variable back within desired limits and 
print out a warning message telling the user that the limit has been exceeded. If optimization is 
not active and a limit has been exceeded, the program will only print out a warning message 
that a limit has been exceeded with the rest of the output for that case. 

1 A VARIABLE SCHEDT IT .TNU 


Many components are configured to automatically use certain types of performance 
maps. The user can input these performance maps in a specified tabular format and the 
program automatically interpolates the tables to obtain the given component performance. A 
good example of this is the compressor and turbine performance maps. The variable scheduling 
component is an extension of this performance map capability. Each variable schedule 

w 8 T r ? mak c e 0n * input value of any of the 11 en ^ ne components (inlets 
through shafts/ to be a function of up to any three other user-specified variables. The value of 

the specification is determined by interpolating a user supplied performance map table. The 
table format is the same as for other component map tables. 
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8.0 CHEMICAL EQUILIBRIUM 


j 1 ? e - de/>ault thermodynamic routine used in the program is preset with properties for air 
and JP4 met (carbon-to-hydrogen ratio of 0.5245 and a heating value of 18600 BTU/lb). For 
many applications, the default thermodynamic routine gives accurate answers If the user 
wants to use a fuel other than JP4 fuel, or operate at high-Mach or high-temperature 
conditions, The default thermodynamic routine may not adequately model these situations, 
■niese possibly significant errors in performance estimates and the dissociation effects are 
discussed in detail in reference 15. 

, ^ chemical dissociation option is included with the program, which enables the program 

to handle just about any combination of chemical species for which thermodynamic data is 
available. The dissociation option uses a rewritten version of the computer program for 
calculation of complex Chemical Equilibrium Compositions (CEC) (ref. 16), to calculate 
thermodynamic properties. The incorporation of the revised CEC program into the NNEP89 
program is discussed further in reference 4. For simplicity, user inputs have been kept to a 
minimum, but still allow a large amount of flexibility for modeling many different types of 
engine cycles. One limitation to the chemical equilibrium option is that the gas stream will 
always be at chemical equilibrium. No chemical kinetics effects will be taken into account, such 
as nozzle gas chemical compositions being frozen at throat conditions. A second limitation to 
the chemical equilibrium option is that, it is only valid in the range of 360 to 9000 °R. 
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9.0 SIMPLIFIED INSTALLATION CALCULATIONS 

The NNEP89 code contains some correlations for estimating installation effects. These 
correlations allow the user to estimate inlet drag or nozzle drag or both. Inlet drags will be 
estimated if variable SPILL is set equal to TRUE. Nozzle boattail drag will be calculated if 
variable BOAT is set equal to TRUE. The program will apply the inlet and nozzle drags to 
calculate the net thrust and total specific fuel consumption (TSFC) with installations drags. 

a ) so P rint t0 the output dataset the calculated inlet and nozzle drags. How to 
use NNEP89 to estimate each drag will be discussed in the next sections. 


The inlet drag, calculated by the code, is the sum of bleed, lip, and spillage drags. 

lgure 6 shows a side view of a 2-dimension inlet, as an example to help define and explain the 
inlet installation drag terms. 



Oblique shock 

Figure 6. Side View of a 2-D Inlet 

Ac denotes the inlet capture area or the total stream tube area for the flow captured or 
affected by the inlet. \ denotes the stream tube area for the airflow entering the engine. A RI n 

area ^ fl ° W captured bv the in,et - but which will be bled or removed from 
the inlet for stability or performance reasons. A LIP is not represented in the figure, but it is the 

stream tube area for the air that impinges on the forward surface of the cowl lip, which will 
cause a drag. A SPL denotes the stream tube area for the flow diverted by the inlet. The 
diversion of this spill flow will cause a drag on the inlet. Each area term can be converted into a 
mass flow by multiplying by the free stream velocity and density and the momentum for each 
term can be determined by multiplying its area by the dynamic pressure, q. 

tntol stJo ^ l aCh ind i vidual ialetdr ?? term , the stream tube area of each term relative to the 
total stream tube area Ac is multiplied by that flow's relative change (or loss) of momentum 

This is called the drag coefficient. C D . The engine airflow, whose stream tube area is denoted by 

A o in f’ i 8 ur « 6, has a ram drag which is always calculated by the program as follows: 


Ram Drag = m * V 0 * 0 031081 


— * ^ tmrnmut itm 


PRECEDING RAM BLANK No, f AKD 



The drag for each other inlet loss is then determined by an equation similar to the following: 

DRAG = C d * Ac*q 


< 


Bleed drag, whose stream tube area is denoted by A BLD , is estimated with two terms, the 

fraction of total capture that will be bled and the change in that flow's momentum. The bleed 
fraction is estimated by the following equation: 

?bld = 0 .016 * (Flight Mach number) 1 - 5 

The change (loss) of momentum is given by the following equation: 

AM bld = (2.0 - (Inlet total pressure recovery) * ^ Flight Mach number ) 

Both equations together give the following relationships for inlet bleed drag coefficient and inlet 
bleed drag: 

Cd bld - ?bld * AM BLD 


and 

DRAGguj = A c * q * C Dbld 

Spillage drag, whose stream tube area is denoted by A SPL> is the total inlet captured airflow 

minus the engine and inlet bleed airflows. The spillage fraction is given by the following 
equations: 


^SPL = = (1-0 - ^BLD ‘ ) 

* 

""TOTAL 

with the limitations: 

if ^SPL > ^SPL = 1*0 
if ^SPL < 00: ^SPL = 0.0 

The change in momentum for the spillage flow is given by one of the following equations: 


if Mach number < 1: AM SPL 


1.33 * (Flight Mach number) 3 ) 
AMINDS 


or 


if Mach number > 1: AM 8PL 



( AMINDS * Flight Mach number )2 


AMINDS is the inlet design Mach number. AMINDS is input by the user, default is 2.7. 




The equations for inlet spillage drag coefficient and inlet spillage drag become: 

Cd spl = *spl * AMspl 
and 

DRAG SPL = Ac * q * C Dgpi 

The cowl lip drag coefficient is estimated using one of the following equations: 

if Mach number < 1: C DLn> = 0.07 - 0.04 * ( Flight Mach number ) 5 


or 

if Mach number > 1: C Dup = 0.07 - 0.04 * ( Flight Mach number - 1.0 )° 25 

If spilling, the cowl lip drag is reduced by multiplying C Dup by ( 1.0 - ? SPL ). Cowl lip drag is 
then calculated using the following equation: 


DRAG LIP = Ac*q*C DLrp 

•n k T ? ta ! ^installation is the sum of bleed, spillage, and cowl lip drags. Inlet drag 
will be calculated if the user sets SPILL equal to TRUE. The user must also inpuUnlet desijf 

thrvariahWr^PT U ? MS S 2 ' 7) ' The USer can input the inlet capture area directly through 
the variable ACAPT or have the program calculate it by setting INLTDS equal to TRUE and & 

inputting a value for the inlet design point spillage fraction in variable SPLDES (default is 0 0) 
TTie program w,U use the present engine parameters, with the inlet design Mach number and 
inlet design spillage fraction to estimate the capture area The user may also modify the 

a , by * he program by in P u tting a value for variable SIZINL. SIZINL is 

teen found to ? * “’ft"!! **** * get 8 new ’ resized capture area. This feature has 

area without ctenHni T determinin e tbe effects of changing the inlet design capture 

area, without changing any other engine or inlet design parameters. 

9.2 NQZZT.F. DR An 

TRITF Th ZZle boattai,dra e is calculated by the program when variable BOAT is set equal to 
™ E , v lere are , no other user inputs for nozzle boattail calculations. Nozzle boattail drae is 
7 U d When nozzle exit area 18 less than the maximum engine nacelle area, as shown in pfgure 



„ a- 

Engine Nacelle^J |a exix 


P^MAXEN 


Figure 7. Engine Drawing for Nozzle Drag Calculations 


35 



i 


T 


T 


Ajjxn. is the nozzle exit area and A maxen i 8 the engine nacelle maximum cross-sectional area 

if^ N arlT u a ? at , CS18 T POint U8ing ** COre fl0W conditions - Th® program assumes that 
if there are multiple nozzles, A^ is the sum of all nozzle exit areas and the nozzles are 

concentric. The boattail drag coefficient, is calculated using the following equations- 


a EXIT 

r = >Wn 


with the limitations: 

*f A r > 1.0: A R = 1.0 

*f A R < 0.05: A r = 0.05 


if Mach number < 1: 


^®boat ^ ) * ( 0.05 + ( Mach number) 16 ) * 0.5 


if Mach number > 1: C D = 0.0476 * — 1,0 * A R 


d BOAT ( Mach number) 1 - 6 

The nozzle boattail drag is given by the following equation: 

drag boat = Amaxen * <J * c Dg 0 A T 


NNEP89 assumes that nozzle boattail drag occurs only if the nn 77 lo „ •*. 

e engine nacelle maximum cross-sectional area. If this is not true if a 


area is less than 


boattail drag is zero, 


*v«b yuan 

is is not true, if Aexjx > Amaxen* nozzle 



► 
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10.0 INPUT/OUTPUT DESCRIPTIONS 


.he^ .„T^rp -ts on unit 9 * * ««° >- 

FORTRAN Namelist input format with &D as thl NailS of the input is read using 

read from Unit 12. Output is written to Unit 10. k name ' Component ma P s are 


IQ .1 GLQBAT. &D NAMTTT j gT INPUTS 




10.1.1 Input/Ontp^ 


AMAC 

DOUTHD 

DRAW 

ITERM 


ITPRT 

LABEL 


LONG 

MAPLOT 

NUMLIN 


-T Write an output data file in the format required by the AMAC 

oZt "ie ^ Sre Writ “" 10 U "“ 14 *" d *PP^X thoZn" 

=F (default) Do not write AMAC data. 

= l ' abe,s ab «ve the DATOUTs on the output 

=* (default) Do not print labels. 

= p r W f ! ) !® ck , diaeram of en £ine on output file. 

=* (default) No diagram will be drawn. 

=0 (default) No terminal printout 

n „l; <«.■ aim.*, 

executing, ' TSFC a " d number of iterations as NNEP89is 

=2^ Same as 1, plus print the number of iterations (NIT) when NIT is a multiple 

is greater^han nS''” ter,m " a ' SCree " Whe " " Urab<,r of lines P ri "’* d 
TsUrZ; N™L P ,N lear termi " al “ ree " Whe " Printed out 

=0 ^ efault [ °o not print tables of component maps. 

-i Flint tables of component maps on output. 

executed. ^ * ' aW “ the top ” f “ P“8o identify the point being 

et^sof WBeZtIuZ^T' ThCn ' in ° be,S 8re deSired fOT ofdesign 

SNoStZwZo^Z ° f lhe C °" VergenCe 11 iS ad ™ M * *» »nve 

-F No printing of convergence history. 

dlr^Zn^raTrr 1 maPS ' includi "« ,he PPetafing points. (See reference 
=F (default) No component map plotting. 

Number of lines available on terminal screen. (Used bv Droeram tr. 

when to automatically clear the terminal „ , ^ p ogram to determine 

naticany clear the terminal screen when using ITERM<0. (Default 
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PINPUT 


PLOT 


TABLES 


XNUM 


value is 20.) Using ITERMcO and NUMLIN is specific to systems that support 
the FORTRAN call CLRSCR. If your system or terminals do not support this 
call, remove the CALL CLRSCR line in subroutine CLRVU. 

=T (default) Print Namelist input for a case on the output prior to the results 
for that case. 

=F Do not print Namelist input . 

=T Enable interactive graphics to draw the engine profile when using the 
WATE89 program. (See reference 17 for more details.) 

=F (default) No WATE89 graphics. 

=T (default) Component maps are used (Unit 12), 

=F Component maps are not used 

Array of characters that will be used to label the flow station on the output. 
Each array location can be up to 4 characters input in Hollerith or quoted string 
format. (Default is the flow station number.) 


ACTL 


ENDIT 


IDONE(i) 


IWAY 


MAXNIT 

MODE 

MODESN 

NCASE 


=0 (default) All controls must be defined and activated by the user. 

=1 Controls to eliminate flow and work errors will be automatically set up and 
will be activated on the first off-design case. Additionally, required controls 
when using CALBLD=T or when using a heat exchanger will be automatically 
activated. Note, because of the many options available, all mixer controls must 
be input by the user. 

=2 Same as 1 above, except the user may interactively vary the independent 
variables chosen to eliminate the flow and work errors. 

=1 Terminate execution after completing the previous case. ENDIT will not 
appear in the global Namelist inputs, but rather in a later input case where the 
user wants the program to terminate. (Default is 0.) 

Array to tell the program if component i has been designed or not. This is used 
to redesipi a single component during the execution of a set of ofT design cases. 
To redesign the entire engine, the user should use IWAY=1 or NCODE=3. 

= 0 Component i is not designed. 

= 1 Component i has been designed. 

Design point switch. (Default value is 1 for the first point, 0 for following points, 
and is internally set to -1 for first pass through an engine at design point.) 
Setting IWAY=1 will execute the next case as a design point, resetting all map 
scale factors. 

= 0 (default) Do not execute WATE89. 

= 1, 2, or 4 Execute WATE89, a weight estimation program. (See reference 10 
for more information.) 

Maximum number of iterations allowed to in order to reach a converged 
solution. (Default is 50.) 

Designates the engine mode the code will use for the next case. (Default value is 

Designates which engine mode is the design mode. (Default value is 1.) 

When set to 1, this tells the program to reset all arrays and read in a new set of 
KONFIG and SPEC information. (In essence, clear everything and get ready to 
start all over again with a new engine configuration and new cycle parameters.) 
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NCODE 


NMODES 

PUNT 


=1 Normal executing. 

=2 Debug mode (print output after each pass) 

PM! 1 through cyc,e ia mada - e * ch »- 

i» rt ' n ' ou ‘ put, and obvia ‘" naad to 

7 4 , Pril \ t p * rtia J derivative matrix each time it is updated and the values of 
independent and dependent variables for each control for each iteration. 

The total number of modes to be configured. (Default value is 1.) 

=T (defauk) Use the values of all SPECs from the last converged case as the 
starting values for the next operating point converged case as the 

-F Use the present values of all SPECs as starting guesses for next point If 

the previous case did not converge, some SPECs may have ridiculous values 

C8Se fr ° m COnverging - ^refore, it is advisable to always 


DEBUG 


NJOPT 


NVOPT 


TOLOPT 


=0 (default) Do not print the debug information when executing optimization 
Print debug information when executing optimization. 

Component number which indicates the location of the dependent variable 
(U d when optimizing. If 0, the dependent variable is not a DATOUT variable.) 
Used when optimizing 

fiSSllSlr ° f 110,7 indiCating WhiCh k the 

" *e°dZ*nde‘ *° * indiCatine Whi ' h DAT0UT ° f NJOPT 

(NOTE: NVOPT >0 minimize dependent variable, 

NVOPT < 0 maximize dependent variable.) 

“St O f oS) ' < ‘ fdePend "' t variable - optimizing, default 


CALBLD 

ELIFE 

NEWEFF 

TMBLAD 
TMVANE 
YE ARB 


too l ing (see fiftfftion 10 . 2.6 nnd reference q % more 

=T Perform turbine bleed flow calculations 
=F (default) No turbine bleed flow calculations. 

m&zt when CALBLD=Tfor ,urbin ‘ biMd «w»«on. 

Actual desired blade metal temperature, ‘R. Default is a function of YEARB. 

Actual desired vane metal temperature, °R. Default is a function of YEARV. 

ear of first service of blade. Only used if TAT ri n-T /» 4 , 

calculations. (Default .due i.7wo” for turbm. bleed 

Yeer of firat service of vane. Only used if CALBLD-T for turhin. hi..n 
calculations. (Default value is 1990.) turbine bleed 


0 . 


10.1.5 Thermodynamic Properties (Including CRCi 


C2HRAT 

FAKRAY(ij) 

ICEC 

OAERAY 

TFUEL 


Carbon to hydrogen atomic ratio, used in duct burner calculations to correct the 
fuel heating value for temperature effects. (Not used if ICEC > 0. Default value 
is 0.5245035801, the value for JIM.) 

Array of fuels to be used by the program in DUCTs and GGENs. j is the 
component number of the DUCT or GGEN. i is an index of reactant number and 
its relative amount to be used as the fuel, up to 6 components. (FARRAY is only 
used if ICEC > 0.) 

= 0 (default) No equilibrium gas properties. 

= 1 Use chemical equilibrium code (CEC) to calculate thermodynamic 
properties. 

Array (i j) of oxidizers to be used by the program in GGENs. j is the component 
number of the GGEN. i is an index of reactant number and its relative amount 
to^be used as the oxidizer, up to 6 components. (OARRAY is only used if ICEC > 

Temperature of fuel added in burner. (Not used if ICEC > 0. Default value is 
530°R.) 


10.1.6 Installation 


ACAPT 

AMINDS 

BOAT 

INLTDS 

INST 

SIZINL 

SPILL 

SPLDES 


Capture area of the inlet, ft 2 - This input is used for the simplified installation 
calculations. This value will be overridden if INLTDS=T. (Default is 0.0). 

Flight Mach number where inlet is designed. Used for simplified installation 
effects calculations when SPILL=T. (Not to be confused with installation effects 
calculated by the INSTAL89 program. Default is 2.7.) 

=T Perform simplified boattail drag calculations. (Not to be confused with 
boattail drag calculated by the INSTAL89 program.) 

=F (default) No boattail drag calculations. 

=T Specifies this case as the sizing point for inlet area. Used for simplified 
installation effects. (Not to be confused with installation effects calculated by 
the INSTAL89 program. Default is FALSE.) 

=0 (default) Do not execute INSTAL89. 

=1 Execute INSTAL89 program. 

Scale factor on the inlet capture area. (Default is 1.0). 

=T Perform simplified installation effects for spillage and lip drag. (Not to be 
confused with spillage or lip drag calculated by the INSTAL89 program.) 

=F (default) No spillage or lip drag calculations. 

Amount of design spillage when INLTDS=T. (Fraction of inlet airflow that is 
spilled to total amount of airflow captured. Default is 0.0). 


This is all of the data that is read on the global Namelist read. Note that comments may 
be included anywhere within the Namelist data. These comments must be preceded by a /* and 
followed by a */. (An example of using this option is given in Section 12.1.) If the ICEC option is 
used, the chemical equilibrium inputs follow the global Namelist read. Refer to Section 10.4 as 
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TART F? T re ^JL f0r C T pl ! te ^ e ? a ?* ° n ? 0W to Mt U P “ in P ut file using the ICEC option. If 
ABLES=T, the code will go to Unit 12 and read in the maps file. See Section 10.3 for a 
description of the format for inputting maps. 


,, J^ r , the c °de reads the global Namelist inputs, it knows how many modes are to be 
read. The code reads m the configuration data (KONFIG array) and specifications (SPEC 
y) or all of these modes. After NMODES number of Namelist data have been read the 
program will execute the MODESN as the design point. For all component types except 
variable control components, the date arrays are read in the following form: 

K0NFIG<1,N)=’NAME',JM1,JM2,JP1,JP2, SPEC(1,N)=V1,V2,V3...V15 

N^E N ™f ™2 m 5?.T"^ U '"^ r .r d **■«"« **• VP. of component, TO, «!»«, for 

^ ’ C ' 8nd the V<dues m *• SPEC arra y are ^fi^d in the following 
sections describing each component. The format is the same as that listed above for variable 

control components except that the SPCNTL array is used at design point to read in the date 
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INPUT 

KONFIGU.N) INLT or 4HINLT or 1 

KONFIGK2.N) Main upstream flow station number (JM1). 

KONFIG{3,N) 0 

KONFIG(4,N) Main downstream flow station number (JP1) 
KONFIG(5,N) 0 


SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 


SPEC(7,N) 


Inlet mass flow, Ib/sec. 

Free stream temperature, °R. (Only used if SPEC(9,N) or ALTP<0.) 

Free stream static pressure, lb /in 2 . (Only used ifSPEC(9,N) or ALTPcO.) 

Inlet additional drag divided by dynamic pressure or table reference number - 
if blank, additional drag is zero Table = F(x,y) where x = Mach number, y = 
inlet exit referred flow, w>/t/p. 

Mach number at inlet. If MACH is input it will override SPEC(5,N). 

Inlet recovery, constant or table reference number. If ETAR is input it will 
override SPEC(6,N). 

=0 Mil Spec is used. 

Mil Spec is: for Mach < 1 Recovery = 1.0 

for Mach > 1 Recovery = 1 - 0.075 (Mach - 1) 1 - 35 . 

Table = F(x,y) where x = inlet exit corrected weight flow, (W\/e/S) y = Mach 
number. 

Maximum permitted flow in table (if SPEC(6,N)= table reference number). 
SPEC(8,N) Scale factor on corrected weight flow = used to read inlet recovery 

table. (If SPEC(6,N)= table reference number.) 

Geometric altitude, feet. If ALTP is input it will override SPEC(9,N). 

Blank 

If nonzero, SPEC(9,N) is geopotential altitude. 

AT to be added to inlet entrance temperature from internal 1962 ARDC 
atmosphere table. Atmosphere table is used if altitude is input. 

Blank 

SPEC(15 d N)°) W ^ ' nlet eXit ’ (ThiS inPUt wil1 override SPEC(l.N) and 
Corrected flow at inlet entrance. (This input will override SPEC(1,N).) 


SPEC(9,N) 

SPEC(IO.N) 

SPECdl.N) 

SPEC(12,N) 

SPEC(13,N) 

SPEC(14,N) 

SPEC(15,N) 


used to read inlet recovery 


OUTPUT 

DATOUT(l,N) Inlet ram drag from table or computed, lb. 
DATOUT(2,N) Free stream velocity, ft/sec. 

DATOUT(3,N) Free stream velocity, knots. 
DATOUT(4,N) Ram temperature ratio. 

DATOUT(5,N) Ram pressure ratio. 


] 
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w 


DAT0UT(6,N) 

DAT0UT(7,N) 

DAT0UT(8,N) 

DAT0UT(9,N) 


Free stream Mach number. 

Intel recovery, ratio of exit total prewar* t„ ram prewure. 

Corrected exit temperature Inlct exit ^wperatun* 

’ 618.67 

Altitude, ft. 


Inlet usage notes: 


respoctted^forL’A Metta^ rcpl,,< * SPECs *• 9 ' "* «• 


for 0 ?T EC(9 ' N) ' to input a value 

pressure, and SPEC(2,N) is ignored ffth^Ii?^V he program 8ets SPEC(3,N)=ambient 
value in SPECe 2 and 3?w !£££ h “ “>!?»• **» P"*™ uw" the 

Md the altitude becomee tew than 0, the IfSPWW “ '•» M«nk, 




10.2.2 Duct or Burner 


INPUT 

KONFIGU.N) 

KONFIG(2 > N) 

KONFIG(3,N) 

KONFIG(4,N) 

KONFIG<5,N) 


'DUCT* or 4HDUCT or 2 

Main upstream flow station number (JM1). 

Secondary upstream flow station («TM2). Used only if duct is used as a passage 
(no burning). 

Main downstream flow station number (JP1). 

Secondary downstream flow station (JP2). Used only if duct is used as a 
passage (no burning). 


SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 

SPEC(7,N) 

SPEC(8,N) 


AP/P, total pressure drop or table reference number. This value is in addition 
to the optional Rayleigh total pressure drop calculated by the program (see 
SPEC(2,N). Table = F(x) where x = burner entrance corrected weight flow 
We/s maD 

W^ S actual ‘ 

Design duct entrance Mach number (optional). This is used at the design point 
to calculate the cross-sectional area of the duct or burner (saved in SPEC(7,N)). 
This area will be used to calculate the Rayleigh total pressure drop from 
heating. 

Additional total pressure drop divided by burner inlet referred flow squared. 

ap/p/(wVt/p) 2 . 

< 0 Desired exit fuel-to-air ratio 
= 0 Duct only 

> 0 Desired exit temperature, °R 

Burner efficiency, ratio of actual heat of combustion to ideal heat of 
combustion, or table reference number (0 if duct). Table = F(x,y) where x = 
burner entrance referred flow, W>/t/p, y = fuel-to-air ratio 

Fuel heating value or table reference number, 0 if not burning. (18300 is 
typical for JP fuel.) - BTU/lb. Table = F(x,y) where x - burner exit 
temperature (SPEC(4,N)), y = burner entrance pressure. 

SeCti ° nal area of du . ct or burner > in2 Calculated at design point if 
SPEC(2,N)>0. This area is used to determine duct or burner inlet Mach 
number. In burners this will also be used to determine Rayleigh pressure 

Ratio of entrance bleed flow to total flow available (duct only). 


SPEC(9,N) 
SPEC(IO.N) 
SPECK 11, N) 


Ratio of exit bleed flow to total flow available (duct only), — 

(*JP2 + fypi) 

Fraction of air not heated (bypassed around combustor portion of burner but 
will be added back before exiting the duct). 

(Only valid if ICEC=1) 

=0 

if SPEC(4,N) > 0 Find fuel-to-air ratio to get desired temperature on fuel-lean 
side 

if SPEC(4,N < 0) Calculate temperature with input fuel-to-air ratio . 

=1 Find stoichiometric fuel-to-air ratio. If entrance fuel-to-air ratio is less 



A 


i 
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I' 


t 


SPEC(12,N) 


SPEC(13,N) 

SPEC(14,N) 

SPEC(15,N) 


than this value, add enough fuel to raise the fuel-to-air ratio to the 
stoichiometric value, otherwise, write error message to output and exit burner. 
=2 Find fuel-to-air ratio to reach desired temperature, but on the fuel-rich 
side. 

=0 No emission calculations. 

=1 Calculate emissions index using a simplified empirical correlation based on 
burner entrance pressure and entrance and exit temperatures. The emissions 
index is given as grams of oxides of nitrogen produced per kilogram of fuel. 

Emissions Index = 0.0070978 * T 4 * ^ *exp p'^j . where P 3 = burner 

entrance pressure, T s = burner entrance temperature, and T 4 = burner exit 
temperature. 

=0 Burner will not reset SPEC(4,N). 

>0 if SPEC(4,N) < Tjm!, set SPEC(4,N) = + 1. 

Blank 

Blank 


OUTPUT 

DATOUT(l,N) 

DATOUT(2,N) 

DATOUT(3,N) 

DATOUT(4,N) 

DATOUT(5,N) 

DATOUT(6,N) 

DATOUT(7,N) 

DATOUT(8,N) 

DATOUT(9,N) 


AP/P, total pressure drop from Rayleigh calculation (if SPEC(2,N) or 
SPEC(7,N) was specified at design point). 

AP/P, total pressure drop specified by SPEC(l.N) + SPEC(3,N). 

Design duct entrance Mach number, equal to SPEC(2,N). 

Ratio of fuel mass flow to duct inlet mass flow. 

Duct cross-sectional area, in 2 . 

Fuel mass flow, lb/hour. 

Duct entrance Mach number (if SPEC(2,N) or SPEC(7,N) is specified at design 
point). 6 

BU I I )n c f5 1C ' enCy (equaI to SPEC < 5 « N ) if ICEC=0, 0 if no burning, 1 if burning 
and ICEC=1). 

Burner outlet temperature before fraction of air not heated (SPEC(10,N) is 
added back to the burner flow, °R. 


Duct or burner usage notes: 


a tr ^ r" USmg the duCt to bypa8S flow » the Program does not account for the pressure 
differences from an entering bypass stream and the main stream. This means that the bypass 

8 pre8sur f 18 only used calculate its enthalpy, total enthalpy will be conserved, but the 
mixed stream will be at the main stream's pressure. 

tb. d n Ct °T. bu "‘ er u compo l nen / t T calculate Mach number and static conditions for 

t e incoming flow if a Mach number (at design-point only) or a cross-sectional area is input. 
Previously, this component would only calculate static conditions if combustion occurred The 
staUc pressure and Mach number for both the main upstream and downstream flows are 
included in the station property output. 

<?Pirm?w E i C(1 nL N) 8h ® u,d be U8ed when a variabl « control varies the burner temperature. 
SPEC(13,N) will keep the burner temperature from dropping below the burner entrance 
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temperature. This situation could cause convergence problems with the program. 
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input 

K0NFIG<1,N) •< 
K0NFICK2.N) B, 
K0NFIG(3,N) 0 
K0NFIG<4,N) N 
K0NFIG(5,N) 0 


'GGEN* or 4HGGEN or 3 

Main upstream flow station number (JM1). 

0 

Main downstream flow station number (JP1). 


SPEC(l.N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 

SPEC(7,N) 

SPEC(8,N)- 

through- 

SPEC(15,N) 


Exit temperature of fuel or fuel and oxidizer stream, °R. 

f“ e ' to oxidizer mass ratio (default = 0.3). Since the fuel-to-oxidizer mass ratio 
is double valued with respect to combustion temperature, an initial guess is 
necessary. 

Gas generator pressure, psia (default is pressure at JM1). 

Fuel mass flow rate, lb/s. 

Oxidant mass flow rate, lb/s. 

Fraction of fuel not included in engine fuel consumption calculations. 

Fraction of oxidizer not included in engine fuel consumption calculations. 


Blank 


OUTPUT 

DATOUT(l,N) 

DATOUT(2,N) 

DATOUTX3.N) 

DATOUT(4,N) 

DATOUT(5,N) 

DATOUT(6,N) 

DATOUTX7.N) 

DATOUTX8.N) 

DATOUT(9,N) 


Gas generator exit temperature, °R. 

Assigned gas generator fuel-to-oxidant mass ratio. 

Gas generator pressure, psia. 

Gas generator fuel mass flow, lb/s. 

Gas generator oxidant mass flow, lh/s. 

n ° W inCluded in fuel 

Calculated fuel to-oxidant mass ratio. 

Total propellant mass flow, Ib/hr. 

ma5S n ° W n<,t ‘ n " Uded “ enSi " e f “ el — »«» 


Gas generator usage notes: 

combustion chamber sending its exhaust into an engine. The user can use it to: d 


U Add a certain amount of a fuel or a fuel and oxidizer stream to 
2) Add a fuel stream at a certain temperature into an engine, or 


an incoming 
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3) Mix specified fuel and oxidizer streams together for a desired temperature or 
mixture ratio (like a rocket iryector and combustion chamber) into the user's 
engine. 

The incoming flow stream determines the gas generator's pressure. If an incomine 
stream does not exist the user must input the pressure. The user specifies in the FARRAY 
aiTay which reactants and their relative amounts the program will use as the fuel. If there m 
also an oxidizer, the user must specify the reactants and their relative amounts in the OARRAY 
array. A discussion and example using a gas generator is included in Section 10.4.2 The gas 

C^ r „re^TtK^ PAIUtAY *° 0Perate; “ 0ARRAY “ ° Wi< ’"* 1 ******* «" 

«. m Gasgenerators are commonly added to a cycle by using a splitter, with the bypass 

“ 5* 6 ga ®f ne ™ tor V*P ut 8 ‘ ream - To model a rocket, flow the user must specify a zero 
bypass on the splitter and exhaust the gas generator flow into a nozzle. 

There are five input combinations depending on the input information the user supplies: 
Input Combination 1: 

ratio SsPErr^n Y ’ ® A1 ^ AY ’ mperature, a starting guess for the fuel-to-oxidant 
ratio (see SP EC(2,N)), and a flow rate for either fuel or oxidizer. The routine will 

calculate the fuel to oxidant mass ratio required to get the desired temperature. 

Input Combination 2: 

U^r suppHes FARRAY OARl^Y, no temperature (or 0), and either a fuel to oxidant 

-rr nd f 7 1 il te I 0r , the I Uel 0r oxidizer; or no fuel to oxidant mass ratio and 
™ tes for t the f “ el and oxidizer. The routine combines the fuel and oxidizer 
m the specified amounts and calculates the exit temperature. 

Input Combination 3: 

O^Y^venT^’," 0 0ARRAY ' ‘ "•* rate far Ih, fad (with no 

^ h , TT” ** y fucl 10 “ r r ‘ ti0 or otidi»r How rate input) 

routano then calculate, the fad', thermodynamic property at the .perilled 
temperature and pressure. 

Input Combination 4: 

User supplies same information as Input Combination 2, except there is incoming flow 
The user must input an oxidizer flow rate if one is desired (fuel to oxidizer mass ratio is 

incoming fT* add# the 8 P ecif,ed amount? of fuel and oxidizer to the 

incoming flow and calculates the mixture s thermodynamic properties. 

Input Combination 5: 

Combination 3, except there is incoming flow. 
The routine then adds the specified amount of fuel, at the specified temperature to the 
incoming flow and calculates the mixture’s thermodynamic properties. 
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input 

KONFIGU.N) *WIN,r or 4HWINJ or 3 

KONFIGK2,N) Main upstream flow station number (JM1) 
KONFIGK3.N) 0 

KONFIGK4.N) Main downstream flow station number (JP1) 
KONFIGK5.N) 0 


SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N)- 

through- 

SPEC(15,N) 


Water-to-air mass flow ratio (maximum value is 0.10). 
Fraction of water vaporized. 

AP/P, total pressure drop. 

Saturation switch: If 0, use SPEC(1,N); if 1, saturate. 


Blank 


O UT PUT 

DATOUT(l,N) 

DATOUT(2,N) 

DATOUT(3,N) 

DATOUT(4,N) 

DATOUT(5,N) 

DATOUIX6.N) 

DATOUT(7,N) 

DATOUT(8,N) 

DATOUT(9,N) 


Actual water-to-air mass flow ratio. 

Input value of fraction of water vaporized. 
Saturation value of water-air mixture. 
Actual fraction of water vaporized. 

AT, change in temperature (°R). 

Water mass flow rate (lb m /hr). 

AP/P , total pressure drop. 

Blank 

Blank 


Water injector usage notes: 

SPFmv^ k" the J Wat f r injector - SPEC(1,N) must be nonzero. The input value of 
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V 


x 


INPUT 

KONFIGO^) 'COMF or 4HCOMP or 4 

KONFIG(2,N ) Main upstream flow station number (JM1). 

KONFIG(3,N) 0 

KONFIGK4,N) Main downstream flow station number (JP1). 
KONFIG(5,N) Bleed flow station for the compressor (JP2). 


SPEC(1,N) R value used to read compressor map tables. 

SPEC(2,N) Ratio of compressor bleed flow to total flow into compressor. 


SPEC(3,N) Scale factor on corrected speed = ' 


(usually input as 1). 


SPEC(4,N) Corrected weight flow (wVe/8) or table reference number. Table = F(x,y,z) 

where x = R value, y = corrected speed (NA/e), z = 3rd dimensional argument 
value. 


SPEC(5,N) Scale factor on corrected weight flow = 


w>/e/8 a 


w — — - ▼rciguL uuw - (usually input as 1). 

SPEC(6,N) Compressor adiabatic efficiency or table reference number. Table = F(x,y,z) 

where x = R value, y = corrected speed (N A/e), z = 3rd dimensional argument 
value. 

SPEC(7,N) Scale factor on compressor efficiency a r| actual map (usually input as 1). 

SPEC(8,N) Compressor pressure ratio or table reference number. Table - F(x,y,z) where 

x = R value, y = corrected speed (NA/©), z = 3rd dimensional argument value. 

PR -1 

SPEC(9,N) Scale factor on pressure ratio = (usually input as 1). 

SPEC(10,N) 3rd dimensional argument value used to read compressor map tables. 

SPEC(11,N) Fractional bleed horsepower lo.’s due to interstage bleed. 

SPEC(12,N) Desired adiabatic efficiency at specified R and N AJe at design point. (Used to 
determine map scale factor for compressor efficiency.) 

SPEC(13,N) Desired pressure ratio at specified R and N/>/e at design point. (Used to 
determine map scale factor for compressor pressure ratio.) 

SPEC(14,N) Corrected speed (NA/e) at design point on maps. 

SPEC(15,N) Blank 


OUTPUT 

DATOUT(l,N) Power required by compressor (negative), hp. 

DATOUT(2,N) Angular velocity of the compressor, rpm. 

DATOUT(3,N) 3rd dimensional argument value used for compressor map tables. 
DATOUT(4,N) R value used for compressor map tables. 

nifliAfmv. m n f Woon/Wcorr..,_ 1 


DATOUT(5,N) Surge margin (percent) * 


PR/PR. 


• 1 x 100. 
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< 


DAT0UT(6,N) 

DAT0UT(7,N) 

DAT0UT(8,N) 

DAT0UT(9,N) 


Corrected speed (NA/e) used for compressor map tables. 
Scale factor on corrected weight flow, (Vfyfo/&). 
Compressor efficiency. 

Compressor pressure ratio. 


Fan or compressor usage notes: 
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10.2.6 Turbine 

input 

KONFIG(l^) 

KONFIG<2,N) 

KONPIG(3,N) 

KONPIG(4,N) 

KONFIGK5.N) 


• -VT 


TURB' or 4HTURB or 6 

Main upstream flow station number (JM1). 

Bleed flow station (JM2) (Matches bleed flow station for the corresponding 
compressor). 

Main downstream flow station number (JP1). 

0 



A 


SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 


SPEC(5,N) 

SPEC(6,N) 


SPEC(7,N) 

SPEC(8,N) 

SPEC(9,N) 


SPEC(10,N) 

SPEC(11,N) 

SPEC(12,N) 

SPEC(13,N) 

SPEC(14,N) 

SPEC(15,N) 


Pressure ratio used to read turbine map tables. 

Ratio of total bleed into turbine to total bleed available. 


Scale factor on corrected speed = 


N />/e, 


actual 


NA/e, 


map 


(usually input as 1). 


Corrected weight flow (ws/e/s) or table reference number. 

Table = F(x,y,z) where x = pressure ratio, y = corrected speed (N A/e), 
z = 3rd dimensional argument value. 


Scale factor on corrected 


weight flow = 


W ^ 8 m«p 


(usually input as 1). 


Turbine adiabatic efficiency or table reference number. Table = F(x t y t z) where 
x = pressure ratio, y = corrected speed (N z = 3rd dimensional argument 
value. 


Scale factor on turbine efficiency 
Scale factor on turbine pressure ratio 


actual^ map (usually input as 1). 

^actual'* 

(usually input as 1). 


Ratio of turbine bleed flow iryected at the turbine entrance to the total bleed 
flow into this turbine. The remainder of the bleed flow into the turbine will be 
added at the turbine exit. 


3rd dimensional argument value used to read turbine map tables. 

Desired adiabatic efficiency at design point. 

NA/e at design point on maps. 

Design point turbine horsepower split (usually set equal to 1). For multiple 
turbines on one shaft, what fraction of the total shaft work required this 
turbine must produce. 

Cooling type code (only used if CALBLD=T, see next page for code definitions). 
Number of turbine stages (only used if CALBLD=T, see next page). 


1 


OUTPUT 

DATOUT(l,N) 

DATOUTX2.N) 

DATOUT(3,N) 

DATOUT(4,N) 

DATOUT(5,N) 


Power produced by the turbine (positive), hp. 

Angular velocity of the turbine, rpm. 

3rd dimensional argument value used for turbine map tables. 
Turbine pressure ratio used for turbine map tables. 

Scale factor on corrected speed, NA/e. 
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, i in mu limp 


DAT0UT(6,N) 

DATOUT^.N) 

DAT0UT(8,N) 

DAT0UT(9,N) 


Corrected speed (nA/5) used for turbine map tables. 
Scale factor on corrected weight flow w>/i/8. 
Turbine efficiency. 

Turbine overall pressure ratio. 


Turbine usage notes: 



The format for the turbine map tables are included in Section 10.3. The program has 
the capability to plot the scaled component maps and the actual operating points on these maps. 
This capability is described in reference 6. 


If performance maps are used, a variable control is required for off design operation of 
the program to drive the flow error on the incoming gas stream to zero. 


X 


I 


I 
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Turbine cooling bleed usage notes: 

CALBLD is set to T when bleed requirements are to be determined. See reference 9 for 
a complete discussion of the turbine cooling analysis. A control must be set to vary SPEC(2,N) 
of the compressor in which the bleed is being removed to drive PERF(15) to zero. Other controls 
may be turned on as required. 



SPEC(14,N) of the turbine is set to indicate the type of cooling row by row. The default 
value is 88 (one stator row and one rotor row with full film cooling). See Table 6 for the 
definitions of cooling types. 


SPEC(15,N) of the turbine is the number of turbine stages and is used only in sizing 
bleed requirements. Default is 1 stage. 

ELIFE = Desired engine life. Default is 10000 hrs. 

NEWEFF = Calculate new turbine efficiency due to cooling. Default is F. 

TMBLAD = Actual desired blade metal temperature, °R. Default is a function of YEARB. 

TMVANE = Actual desired vane metal temperature, °R. Default is a function of YEARV. 

> YEARB = Rotor blade first year of service. Used to determine blade temperature Default is 

l , 1990. 

i YEARV = Stator vane first year of service. Used to determine vane temperature. Default is 

! 1990. 



Table 6. Cooling Type Definitions. 
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INPUT 

KONFIGU.N) 'HTEX’or 4HHTEX or 6 

KONFIG(2 > N) Upstream How station number for the flow stream being heated (JM1). 
KONFIGK3.N) Upstream flow station number for the flow stream being cooled (JM2). 
KONFIG(4,N) Downstream flow station number for the flow stream being heated (JP1). 
KONFIG(5,N) Downstream flow station number for the flow stream being cooled (JP2). 


SPEC(1 J N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 


SPEC(5,N) 


AP/P or table reference number for heated flow stream (JM1). 

Table = F(x) where x = Referred flow of JM1 stream, W\/t/p. 

AP/P or table reference number for cooled flow stream (JM2). 

Table = F(x) where x = Referred flow of JM2 stream, W\/t/p. 

Guess value for temperature rise of main stream (JM1), °R. 

Heat exchanger effectiveness (ratio of actual thermal energy transfer to ideal 
thermal energy transfer) or table reference number. 

Table = F(x,y) where x = weight flow of the stream being heated, W^, y = 
weight flow ratio of the stream being heated to the stream being cooled, 

WjMl 

Scale factor on heat exchanger effectiveness. 


SPECs 6 through 13 are only active if ICEC=1 (see Heat exchanger usage notes). 
SPEC(6,N) Specific enthalpy for heated flow stream (JM1), BTU/lb. 

SPEC(7,N) Ratio of specific heats (gamma) for heated flow stream (JM1). 

SPEC(8,N) Gas constant (R) for heated flow stream (JM1), BTU/lb-°R. 

SPEC(9,N) Specific enthalpy for cooled flow stream (JM2), BTU/lb. 

SPEC(IO.N) Ratio of specific heats (gamma) for cooled flow stream (JM2). 
SPECQ 1,N) Gas constant (R) for cooled flow stream (JM2), BTU/lb-°R. 
SPEC(12,N) Temperature for heated flow stream (JM1), °R. 

SPEC(13,N) Temperature for cooled flow stream (JM2) °R 
SPEC(14,N) Blank 

SPEC(15,N) Blank 


OUTPUT 

DATOUT(l,N) AP/P, total pressure loss for heated flow stream (JM1), equal to SPEC(l.N). 

DATOUT(2,N) AP/P, total pressure loss for cooled flow stream (JM2), equal to SPEC(2 N) 
DATOUT(3,N) Blank ’ ' 

DATOUT(4,N) Heat exchanger effectiveness, equal to SPEC(4,N). 

DATOUT(5,N) Scale factor on effectiveness, equal to SPEC(5,N). 

DATOUT(6,N) Calculated temperature rise of heated flow stream (JM1),®R. 
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DAT0UT(7,N) 


DAT0UT(8,N) 

DAT0UT(9,N) 


Temperature rise of heated stream divided by the difference between the 
entrance temperatures of the heated and cooled streams 

(T JM1 T JM2> ' 

Temperature rise difference, ((guess value/calculated value)-l). 
Calculated temperature rise of cooled stream (JM2),°R. 


Heat exchanger usage notes: 

exchanw! Sw P ° int ’ ** a contro1 to Vdl 7 either the estimate for the heat 

changer effectiveness or the temperature rise of the main stream such that these two uppr*. 

SK 1 W » — to "design" the heat Tha„£ r ^SSSj 1 S°* 

£ ™ tr ° l8 ar V Urn ^ d J ° ff ’ S ° 88 n0t to change ofcher eye 1 ® Parameters whfie designing the 
heat exchanger. For off-design points, this control is turned off. resigning the 

• ... Tb heat f. d and coo,ed streams can be reversed without reconfiguring the engine bv 

inputting a negative temperature rise (SPEC(3,N) < 0 ). ^ g ‘ gme by 

* PEC ° 6 through have bee* 1 added to override the calling of the CEC routines to get 
thermodynamic properties for one or both of the upstream flows (either JM1 or TM2 

The CFC ? ^ allows the user to input the thermodynamic properties for that stream 

The CEC calls are overridden for the JM1 (main stream) if SPEC(7 N) > 0 The CEC rails 
overridden for the JM2 (secondary stream) if SPEC(10,N) >0. 



10.2.8 Flow Splitter 


MUZ 

KONFIG(l,N) ’SPLT or 4HSPLT or 7 

KONFIG<2,N) Main upstream flow station number (JM1). 

KONFIG(3,N) 0 

KONFIG(4,N) Main downstream flow station number of the split flow (JP1). 
KONFIGK5,N) Secondary downstream flow station number of the split flow (JP2). 


SPEC(1,N) 


SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N)- 

through- 

SPEC(15,N) 


Bypass ratio, ratio of secondary downstream mass flow to main downstream 

- A JP2 
mass flow, . 

*jpi 

AP/P of the main flow stream (JP1). 

AP/P of the secondary (bypass) flow stream (JP2). 


Blank 


OUTPUT 

DATOUT(l,N) 


DATOUT(2,N) 

DAT0U1X3,N) 

DATOUTX4.N) 

DATOUT(5,N) 

DATOUT(6,N)- 

through- 

DATOUT(9,N) 


Bypass ratio, ratio of secondary downstream mass flow to main downstream 

mass flow, . 

*JPi 

AP/P in the main flow stream (JP1). 

AP/P in the secondary (bypass) flow stream (JP2). 

Total weight flow of primary stream (JM1), Ib/s. 

Sum of weight flows of reactants of primary stream (JM1), lb/s. 


Blank 


Flow splitter usage notes: 


Each splitter bypass stream must eventually end in a nozzle or mixer. If such is not the 
case, use a DUCT instead of a splitter, setting SPECO.N) of the DUCT to the desired bypass 
ratio. The program also follows the secondary (bypass) stream through other components (if 
any) and keeps track of the bypass stream's pressure (as opposed to how the program handles 
bypass streams from DUCTs ). This is discussed in Section 5.1.2. 
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EF 

10.2.9 Miiet-fli 


TciUlF 


KONFIGU.N) 
KONFIG(2,N) 

XONFIG<3,N) 

KONFIG(4,N) 
KONFIG(5,N) 

SPEC(1,N) 
SPEC(2,N) 
SPEC(3,N) 

SPEC(4,N) 
SPEC(5,N) 

SPECK6.N) 
SPEC(7,N) 

SPEC(8,N) 
SPEC(9,N) 

SPECX10,N) 
SPECX11.N) 

SPECK 12, N) 

SPECK 13, N) 
SPECK 14^0 

15 


* — — — 


Eiector 



■MIXR' or 4HMIXR or 8 

Main upstream flow station number (JM1) if SPEC(7,N) = 0. 1 

Ejector secondary flow stream if SPEC(7,N) > 0. 

Secondary upstream flow station number (JM2) if SPEC(7,N) = 0. M 

Ejector primary flow stream if SPEC(7,N) > 0. j 

Main downstream flow station number (JP1). 

0 


Mixer inlet area of JM1 stream, in 2 . Calculated at design point using 
SPEC(3,N) and design point operating conditions. 

Mixer inlet area of JM2 stream, in 2 . Calculated at design point using 
SPEC(3,N) and design point operating conditions. 

(Required at design point to determine mixer inlet areas.) 

< 1 Mach number of JM1 stream at design point 
> 1 Total to static pressure ratio of JM1 stream at design point 

If SPEC(4,N) >0 momentum coefficient for mixer exit stream (JP1); ignore 
SPECU3.N) and SPEC(14,N). 

=0 Independent mixer inlet areas. 

=1 Total mixer inlet area is held constant, any change in the secondary (JM2) 
area has a corresponding change in the primary (JM1) area. 

=0 Inactive. 

=1 (Active at design point only). SPEC(3,N) will be the minimum total to 
static pressure ratio or minimum Mach number of the stream that has the 
lowest total pressure (either JM1 or JM2). 


=0 Mixer will mix two subsonic streams. 

=1 Supersonic mixer. At design point, expand the JM2 stream to the same 
static pressure as the JM1 stream. Off-design, do the same calculation as 
SPEC(7,N) = 2. 

=2 Do ejector calculation (JM2 stream is supersonic). (See mixer usage 
notes.) 


= 0 Subsonic solution to mixing equations. 

= 1 Supersonic solution to mixing equations. 


If SPEC(7,N) > 0, design point ejector primary (JM2 stream) Mach number 
(typically around 1.5 to 3), otherwise blank. 


Mixer area ratio = 


Mixer exit mw 

Tout mixer inlet wee’ 


Injection angle (degrees) for 
JM1 if SPEC(7,N) = 0, 
JM2 if SPEC(7,N) > 0. 
(Default is 0.0) 


Factor used to check for over-expansion in ejector stream (JM2) if SPEC(7,N) > 
0. (See Mixer usage notes. Default is 0.33). 

If SPECK4.N) S 0 Momentum coefficient for JM1 stream, otherwise blank. 

If SPECK4,N) s 0 Momentum coefficient for JM2 stream, otherwise blank. 


0 
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SPEC(15,N) Leave blank. 


OUTPUT 

DATOUTd,N) Main stream flow area, (JM1 stream), in 2 . 

DATOUT(2,N) Secondary stream flow area, (JM2 stream), in 2 . 

DATOUT(3,N) Total to static pressure ratio of JM1 stream. 

DATOUT(4,N) Total to static pressure ratio of JM2 stream. 

DATOUT(5,N) Velocity of JM1 stream, ft/s. 

DATOUT(6,N) Velocity of JM2 stream, ft/s. 

DATOUT( 7 , N) Exit mixed flow velocity, ft/s. 

DATOUT(8,N) Static pressure difference between the JM1 and JM2 streams divided by the 
average static pressure. 

DATOUTX9,N) Exit area/total inlet area. 

Mixer usage notes: 

A new, improved mixer routine has been put into NNEP89. The current mixer solves 
the equations of mass, momentum, and energy. The current mixer will increase the exit area, if 
necessary, to have a valid solution to the mixing equations. Mathematically, the mixer inlet ’ 
streams do not have to be at the same static pressure to have a valid mixed solution unlike the 
original NNEP mixer subroutine. The original mixer subroutine assumed a constant static 
pressure for all mixer entrance and exit streams, unless the ejector option was used. 

The original velocity coefficient term has been replaced with more realistic momentum 
coefficient terms. These terms are used in the conservation of momentum equation. A velocity 
coefficient term is only applied to a flow stream's velocity contribution to momentum. The 
current momentum terms are applied to a flow stream's velocity and pressure contribution. 

The current mi/er also includes the exit static pressure and Mach number in the station 
property printout for all gas streams. If the current mixer has to increase the mixer exit area to 
get a valid solution to the mixing equations, a warning message will be included in the output. 

If SPEC(? ,N) > 0, the mixer will do ejector calculations. The ejector primary (JM2) 
throat area will be determined from the design point conditions. The JM2 stream's Mach 
number is input at design point and determines its exit-to-throat area ratio (like a C-D nozzle) 
Off design, the throat area and area ratio, with its inlet total conditions will determine the 
ejector primary's mass flow and inlet static conditions. Any time the ejector primary static 
pressure (JM2) is less than SPEC(12,N)*ejector secondary static pressure (JM1), the mixer 
subroutine assumes that the ejector primary flow (JM2) separates. A message warning of this 
condition will be included in the output. The subroutine then recalculates the ejector primary 
stream's static conditions such that its static pressure is equal to the ejector secondary stream's 
(JM1) static pressure. The ejector primary velocity and flow area consistent with its new static 
pressure will be used for the remainder of that ejector calculation. 
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10.2.10 uau i» 


TcTww 


KONFICK1W NOZZ or 4HNOZZ or 9 

KONFIG<2 r N) Main upstream flow station number (JM1) 
KONFIG<3 t N) 0 

KONFIG(4,N) Main downstream flow station number (JP1) 
KONFKXSjN) 0 


SPEC(1.N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 


SPEC(7,N) 

SPEC(8,N) 

SPEC(9,N) 

SPEC(10,N) 


SPEC(11,N) 


“f 1 f ° r C ° nver p nt ’ J thr0at for CD The flow area is 

calculated at the design point based on present cycle conditions. 

Flow or discharge coefficient (C D ) or table reference number. 

Table = F(x) where x = nozzle overall pressure ratio. 

Nozzle exit static pressure, psia. (Default is ambient pressure). 

Ambient static pressure, psia. If 0, see SPEC(9,N). 

Velocity coefficient (C v ) or table reference number. Table = F(x,y) 

where X = c ” tr * nc< 101,1 _Noarte eat sre. 

Noule throat static pressure » ? = Nozzle throat area • 

Geometry Switch: 

=0 Convergent nozzle. 

wseJTsPECn^^n^^ °jP* n8ion to ambient 8tatic Pressure for all 
lf „ > °> then perfect expansion to SPEC(3,N). 

ThateaT i a a n0ZZ *” ^ ozjde e ” 4 throat area ratio is calculated at design 
Nozzle variable area switch: 

-1 t °* the in ? Ut ^ Ue 8pecified b y SPEC(1,N) at off-design. 

. , 7 ^ throat area for C-D nozzles or exit area for convergent nozzles to 
match the area required to pass all the flow entering the nozzle* 

Nozzle exit to throat area rrtio. Used at design point and off-design for Tvne 3 
nozzles. Also caloilated for Types 1 and 2 noTzles at design poSfSut on* 

Sr^e itoilTs U “ 41118 ‘ ?r ° ff deWgn nozzle calcuJa4i °ns. It will be ignored 

SPFm 40 com P° n ent number of the inlet. If 

SPEC 3 UfAefrL ? ’ a " ! e “ tha " ° r equal 40 °’ Ae P ro 8 ram will use 
staUc pnLw* pre88ure of inI ® 4 1) ^ the nozzle ambient 

C-D (nozzle Types 1, 2, or 3) exit divergence half angle for divergence loss 

* ng J e , may *" ** exit nngufof the ,hro?d, theplug(if 
y), or both the shroud and plug, depending on the nozzle type switch for 

tS3SL u« 

length for diverge. In. calculntiona, in. Whan 
**»*•"« half ntf. (SPEC(IOJO) la calculatad from tha km«n 

in SPECnn*Kn *nt**i, 2 1 '* y* 1 * * al uk * P™»danca over my mf le cpaciflad 
in aPEC(lO.N). Blank if n« divergence loaa calculation ia deeired. 
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SPEC(12,N) 


SPEC(13,N) 

SPEC(14,N) 

SPEC(15,N) 


5TdSn7i e K^- 8WitCh f0r C ; D nozzle loss calculations: 

-(Mdefault) disables divergence calculations, divergence loss coefficient set to 

=1 Axisymmetric C-D nozzle. 

=2 2-D exit C-D nozzle. 

tro^ aymmetTiC PlUg C '° n ° Mle with equal wal1 ejrit half angles on plug and 

shroud^ CXit PlUg °’ D n0 “ le ™ th 6qual wal1 exit haJf angles on plug and 

-6 C x'2 ^°“ le 7^ cylindrical (non-divergent) shroud, 

wtllf. D 1 PUe (Wed8C) C *° n ° zzle ^ P ar aUel (non-diveigent) shroud 

Blank 

Blank 

Blank 


OUTPUT 

DATOUT(l,N) Gross jet thrust, lb. 

DATOUT(2,N) Actual jet velocity, ft/s. 

DATOUT(3,N) Total to static pressure ratio at throat. 

DATOUT(4,N) Nozzle exit area, in^. 

DATOUT(5,N) Nozzle throat area, in*. 

DATOUT(6,N) Flow or discharge coefficient, C D . 

DATOUT(7,N) Velocity coefficient, Cy. 

DATOUT(8,N) Critical pressure ratio at throat. 

DATOUT(9,N) Overall nozzle pressure ratio: nozzle entrance total to exit static. 

Nozzle usage notes: 

°! NNEP ri for *-™&<**™*i>* 

thrust or a pressure drag term nrnnf h* innh a a • ^ am ^ ient pressures, a pressure 

SSSsTJSC £S2~“ "“ I,w h '" ce ’ 

that losses due to expansion to other than ambient rnr»d*t° VereXpan ** e< * ° r undere *P and ®d- Note 
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coefficients and/or discharge coefficients. These loss coefficients are only applied to the 
momentum thrust calculations. 

The fixed area ratio calculations are, invoked only when certain inputs are specified. 

Thus input files written for earlier versions of NNEP will produce the same results as before. 

Hie two methods for fixing the nozzle exit to throat area ratio are: 

1. Set SPEC(6,N) = 2 (a "Type 2" C-D nozzle). For Type 2 C-D nozzles, the nozzle exit to 
throat area ratio required for perfect expansion is calculated at the design point 
That area ratio is stored in SPEC(8,N) and is used for all subsequent off-design cases. 
Unless, of course, the user specifies a control to vary the area ratio or inputs a new 
area ratio in SPEC(8,N). 

Note: the user, while exercising caution, may at any time "redesign" a Type 2 C-D 
nozzle by setting IDONE(N) = 0 or by setting SPEC(6,N) = 1 for that case and then 
resetting SPEC(6,N) = 2. 

2. Set SPEC(6,N) = 3 (a "Type 3" C-D nozzle). For Type 3 C-D nozzles, the user specifies 
the desired area ratio by supplying a value in or putting a control on DATINP(8,N). 
The area ratio in DATINP(8,N) is used exclusively in all nozzle thrust calculations. 


Another addition to the nozzle routine is the ability to apply divergence fi.e., non-axial) 
gross thrust losses to converging-diverging nozzles. This is discussed further in reference 6. 
These losses are accounted for by applying an analytically derived, geometry dependent 
divergence loss coefficient to the gross momentum thrust. The user is offered a choice of six 
different C-D nozzle geometries (shown in Table 7) for divergence loss calculations. 








+ 
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10.2.11 Load or PronelW 


INPTJT 

KONFIG(l,N) 

KONFIGK2.N) 

KONFIGK3.N) 

KONFIGK4.N) 

KONFIG5.N) 


LOAD' or 4HLOAD or 10 
0 
0 
0 
0 


SPEC(1,N) 

SPEC(2,N) 


SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 

SPEC(7,N) 

SPEC(8,N) 


SPEC(9,N) 


SPEC(IO.N)- 

through- 

SPEC(15,N) 


Horsepower (negative for load or positive for power input to cycle) or table 
reference number. Table = F(x) where x = actual shaft rotational speed, rpm. 

Propeller efficiency ( shaft power input or table reference number. Note: 
loads other than propellers may be connected to shafts; for such loads only 
SPEC(1,N) is specified, all others must be blank. Table = F(x,y,z) where x = 


coefficient of power, , y 

Mach number. 


= advance ratio, 


It • Flight velocity 
Tip speed 


z = free stream 


Propeller thrust divided by shaft horsepower for static thrust calculations 
when tables are not used, lb/HP. 

Design point power loading (shaft horsepower divided by propeller diameter 
squared). 

Design tip speed of the propeller, ft/s. 

Blank 


Blank 


A constant scaling factor which is multiplied by the results taken from the 
propeller map. 'Hus input scales the efficiency, thrust, and coefficient of thrust 
when propeller maps are specified in SPEC(2,N). SPEC(8,N) will override 
SPEC(9,N) if both are input non-zero. 

Desired propeller design point efficiency when using propeller maps. This 
input determines a scale factor which will be applied to all subsequent off- 
design cases. SPEC(8,N) will override SPEC(9,N) if both are input non-zero. 


Blank 


QUTPUT 

DATOUT(l,N) Load power, hp. 

DATOUT(2,N) Actual shaft rotational speed, rpm. 


DATOUT(3,N) Propeller thrust ^ lb - <1™^ velocity is zero, the equation 

for propeller thrust becomes indeterminate and will be set to zero if SPEC(3 N) 
is not input or tables are not used.) 


DATOUT(4,N) 


Advance ratio, 


K • Flight velocity 
Tip ipeed 
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DAT0UT(5,N) Coefficient of power, 

p n* D 5 ’ 

DATOUT(6,N) Coefficient of thrust,- 2 ^ 

p n* D 4 ' 

DATOUT(7,N) Propeller tip speed, ft/s. 

DATOUT(8,N) Propeller diameter, ft 

DATOUT(9,N) Adjusted propeller efficiency (using SPEC(8,N), or SPEC(9,N)). 
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10.2.12 Shaft 


icraw 


KONFIG{l,N) 

KONFIG(2,N) 

KONFKX3.N) 

KONFIG(4,N) 

KONFIGK5.N) 

SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 

SPEC(6,N) 

SPEC(7,N) 

SPEC(8,N) 

SPEC(9,N) 

SPEC(10,N)- 

through- 

SPEC(15,N) 

output 

DATOUTd.N) 

DATOUT(2,N) 

DATOUT(3,N) 

DATOUT(4,N) 

DATOUT(5,N) 

DATOUT(6,N) 

DATOUT(7,N) 

DATOUT(8,N) 

DATOUTO.N) 


SHFT or 4HSHFT or 11 

Component number (JM1) of the first component coupled to shaft. 

Component number (JM2) of the second component coupled to this shaft. 
Component number (JP1) of the third component coupled to this shaft. 
Component number (JP2) of the fourth component coupled to thin shaft 

Actual shaft rotational speed, rpm. 

Gear ratio of JM1 component, component rp m/shaft rpm. 

Gear ratio of JM2 component, component rpm/shaft rpm. 

Gear ratio of JP1 component, component rpm/shaft rpm. 

Gear ratio of JP2 component, component rpm/shaft rpm. 

Mechanical efficiency of JM1 component, actual power/ideal power. 
Mechanical efficiency of JM2 component, actual power/ideal power. 
Mechanical efficiency of JP1 component, actual power/ideal power. 
Mechanical efficiency of JP2 component, actual power/ideal power. 


Blank 


Ihafl^hp P °* er ' ‘ he SUm ° rDATOUTtl ' N ' ° f al > components connected to the 


Actual shaft rotational speed, rpm. 

Actual rotational speed of component JM1, rpm. 

Actual rotational speed of component JM2, rpm. 

Actual rotational speed of component JP1, rpm. 

Actual rotational speed of component JP2 rpm 
Blank 

Ratio of net shaft power to total power. Total power is the 
values of all loads on this shaft divided by two. 

Blank 


sum of the absolute 


Shaft usage notes: 


work wT DATOim'S’i.? 1 ?" Uire<i » M ° P *I'* ti0 " of “ drive “>« .hrfl 

, ' U ATOUT(8,N) ) to zero. If one shaft is to be connected to another shaft in nrH»r 

fSwTnT 5 t ] ian ^5°7P° nenU on »m« shaft, the user must configure the shafts using the 

— ° flh * — 
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1Q.2.13 Variahlft 


to the program airflow otw wk^or^ the U8er define as in Put 

specify that some components operate at a a *?* ^ nve ^° 2 ®ro. The user may also 

at a constant surge margin). Variable controlei^cifr to “mpreMora 
V ‘"* bl '* “ ' l,an8e to «>» ■»<*.! te iSSto in £e 




KONFIGK1.N) CNTL' or 4HCNTL or 12 

KONFIG{2,N) 0 

KONFIGK3.N) 0 

KONFIGK4,N) 0 

KONFIG{5,N) 0 


is used for input at «£dl^ The SPCNTL variable 

ensign ease, SPCNTL values do not re tain 

The first two inputs describe the independent variable. 

SPCNTL^^N) tJT SPEC nUmber ° f the independent variable. 

Thenl^^ 

PCNTL(3,N) Type of dependent variable being specified. 

-0 if the depend^^Tf V ‘ ri ‘ b, « tf “ “ * DATO “- 

SWNn.g.N, Value of to. dependent 

^ ' *• *** «*»■ 

= 1 Default value of 0.001 is used (0.0001 if NVOPT * 0) 

Minimum allowable value of independent variable - (if 0 f hi. • * • • 

Maximum allowable value of independent variable - (if 0, this 

®5to«SS^ h* U8ed in P^ce of 

I ^n°L variabl ® h®* 11 * «P«cified as a target value 

= I DATOuT^' 


SPCNTU7.N) 

SPCNTU8.N) 

SPCNTL(9,N) 

The next three 
SPCNTU10.N) 



1 
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SPCNTL(11,N) 

SPCNTL(12,N) 


= DINP* or 4HDINP , or 300 - if the target value is a SPEC. 

= PERF or 4HPERF, or 400 - if the target value is a performance property. 

Number that specifies the station property, DATOUT, SPEC or performance 
property of SPCNTLdO.N). 

= Flow station number if target value is a station property. 

= Component number if the target value is a DATOUT or SPEC. 

= 0 if the target value is a performance property. 


SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N) 


SPEC(6,N) 


OFF-DESIGN SPff.On 

b°de^ed asfonows° ntr °^ Specification8 can onl y 1)6 modified using the SPEC array. This array 

SPEC(1,N) Fraction of target value used for marching (see Marching explanation below). 
Minimum allowable value. See variable control usage notes. 

Maximum allowable value. See variable control usage notes. 

The SPEC number of the independent variable. 

Value of the dependent variable to be achieved, 

(Only if SPCNTL( 10, N) at design point or SPEC(IO.N) at off-design equals 0). 
Station property number if dependent variable is a station property, otherwise 

DOUT number if dependent variable is a DATOUT, otherwise blank. 

Performance property number if dependent variable is a 'PERF, otherwise 
blank. 

Tolerance on value to be achieved. If zero, control is off. If SPEC(9,N) = 1 
default value of 0.001 is used (0.0001 if NVOPT * 0). 

SPEC(10,N) thro “ghSPEC(12,N) have the same value as SPCNTL(10,N) through 
orCNTIX 12, N) 

Type of dependent variable, (STAP, DOUT, DINP or PERF). 

Number that specifies the station property, DATOUT, SPEC or performance 
property for a variable target value. 

Flow station number or component number for the variable target value. 

Blank 

Blank 
Blank 


SPEC(7,N) 

SPEC(8,N) 


SPEC(9,N) 


SPEC(10,N) 

SPEC(ll.N) 


SPEC(12,N) 

SPEC(13,N) 

SPEC(14,N) 

SPEC(15,N) 


OUTPUT 

There are no DATOUTS or output for CNTLs. 


Variable control usage notes: 


v.lu.. S" variable., the input maaimum and minimum 

a es will not have any aiTect on the value of the independent variable uaed bv the nrnnam TO 

Uw uMr input ’ ““ ~ 

essage in the output file. If the user is using optimization and the program exceeds a limit 
the program will enalize the optimization function for exceeding the limits That penalty mav 
prevent the program from exceeding the limits. P y may 




I 
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The array of available station properties are: 


Definition 

Weight flow, Ih/s 
Total pressure, psia 
Total temperature, 

Fuel-to-air ratio 

Referred Aow,W\/t/p Note: corrected flow, w>/e/8, is 
printed on the output W\/t/p = 1.5497 w>/e/8 
Mach number 
Static pressure, psia 

. , _ , _ Wcorr, - WcoiTo 

Interface corrected flow error = ■j — 

2 (WcoiTj + Wcorr 2 ) 

where Wcon^ is the corrected flow out of the component 
immediately upstream of the flow station and 

Wcorr 2 is the corrected flow out of the component immediately 
downstream of the flow station. 


The array of available performance properties are: 
Performance Property Nnmhpy Definition 


Total engine airflow (Ib/s) 

Gross jet thrust (lb) 

Fuel flow (lb/hr) (includes oxidizer flow from gas generator) 

Net jet thrust (lb) 

T~rr Fuel flow (lh/hr) , 

~ Net jet thrust (lb) ’ 

Net thrust/airflow , lb-s/lb 

Total inlet drag, lb 

Total brake shaft horsepower 

Net thrust with installation drag, lb 

Net TSFC with installation drag, lb/lb-hr 

Inlet drag flip, bleed and spillage), lb 

Boattail drag, lb 

Blank 

Blank 

Sum of SPEC(2,N) of all turbines - one. See turbine cooling bleed 
calculations for usage. 

Emission Index, grams NO x per kilogram fuel 
Total shaft horsepower of all propellers 
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18 not equjJ to zero > then a variable target value wUl be used instead of 
SPCNTU6.N). SPCNTL(10,N), SPCNTU11.N), and SPCNTL(12,N) are used to allow the user 
to have the program drive the dependent variable equal to the value or the multiple of a 
specified station property, data output, data input (SPEC), or performance property. 

1116 following example demonstrates one use of variable target values. Suppose one has 
a turbofan with 2 nozzles and would like them to have equal jet velocities. The main nozzle is 
component 8, the secondary nozzle is component 20. Just before the secondary nozzle is duct 19 
much is used as an afterburner. One can now vary the temperature (SPEC 4) of the 
afterburner (DUCT 19) to get the jet velocity (DATOUT 2) of nozzle 20 equal to the jet velocity 
ot nozzle 8. This must be done in this manner because we do not know "a priori" the jet velocity 
of nozzle 8. For this example, the control will be number 41 and will have a tolerance of 0.002. 
The minimum and maximum values for the burner temperature are 800 and 3600 °R, 
respectively. The control would be input as follows: 

KONFIGd > 41)='CNTL’,SPCNTL(l,41)=4,19, , DOUT,2,20,0,0.002,800,3600,4HIX)UT,2,8, 


Marching is a special option on variable controls which allows the user to reset the 
target value ( SPEC(5£T) ) of the variable control to a multiple of that valuXm Se pilous, 
converged case. Marching may be initiated in any off-design case by setting SPEC(1 N) of the 
appropriate variable control to the desired multiplying factor that is to be useJ for the 
marching. As long as SPEC(1,N) of the variable control is not zero, the multiplying factor will 
be applied m all subsequent cases. In each case the value of the dependent variable from the 
previous converged case is multiplied by the SPEC(1,N) value and this new value is used as the 
target value of the variable control. 

In If a variable target value is used, marching will work differently than described above 
In toiscasetheyanable specified as the target (SPCNTL(10,N) through SPCNTU12 N)) is 
multiplied by SPiiCd.N) As the program iterates to a converged solution, this result is 
continuously updated to become the new target value of the variable control. 

„ n , . fo [ lowing sample demonstrates the use of marching. Suppose one wishes to make 

T\ T? 0 / Ver8US n0zzle area at Mach 14 and 40000 feet. A variable control 
component is included defining nozzle area as the independent variable and thrust to airflow 

. .? as the depe ™}e nfc variable. An actual target value is not input because it will be computed 

numil f^f m ™ ,8C0 " tr0 J 18 tu ™ d off for the design case. Assuming the component 
number of the nozzle is 10, the variable control would be specified as follows: 

KONFIGK 1,20)='CNTL', SPCNTU 1,20)= 1, 10,'PERF,6,0,0,0 

N ow the engine design point is executed as well as any desired off-design cases. At the point 
where the user wishes to begin marching, the variable control is turned on and the marching 
value is set to its desired value (1.1 in this example) as follows: marching 

AD SPECU,20)=1.1,SPEC(9,20)=1, 4 END 


0 . 
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. . . ™ ™ a ^. fol,owed by as many subsequent cases as desired. When the program 
detects SPECO20) is not equal to zero, it will take the last value of the dependent variable (in 
this case PERF(6)), multiply it by SPEC(1,N), and store that value in SPEC(5) of the variable 
control. The control will then drive the dependent variable to this new target value. 

Marching can also be used when there is a variable target value. Using the previous 
M.mpte .flh* variable target value from Section 10.2.13.1, suppose that the uLr wants the jet 

« ZI J k* 135 times the ve,ocit !' of the Primary nozzle. Ihe user 
would set SPEC(1,41)=1.35 and the program would vary the burner temperature of duct 19 to 
|8t the jet velocity (DATOUT 2) of nozzle 20 to 1.35 times the jet velocity (DATOUT 2) of nozzle 
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10^.14 Variable Ontimiratinn 


IcIuIF 


KONFKKl^) 

KONFIG(2,N) 

KONFIG(3,N) 

KONFIG(4,N) 

KONFICK5.N) 


'OPTV or 4HOPTV or 13 
0 
0 


Component number containing the independent variable. 
0 


SPEC(1,N) 

SPEC(2,N) 

SPEC(3,N) 

SPEC(4,N) 

SPEC(5,N)- 

through- 

SPEC(8,N) 

SPEC(9,N) 

SPEC(10,N)- 

through- 

SPEC(15,N) 


Blank 

Minimum allowable value of the optimized dependent variable (if = 0, there is 
no minimum constraint). 

Maximum allowable value of the optimized dependent variable (if = 0, there is 
no maximum constraint). 

A value from 1 to 15 indicating which SPEC of component KONFIG(4,N) is the 
independent variable. 


Blank 

Switch to turn the OPTV variable on or off. 
If set = 0, the variable is off. 

If set = 1, the variable is on. 


Blank 





. There are additional inputs to NNEP89 when OPTV components are present OPTVs 
te the program which variables are the independent variables. The next variables tell the 
program whether or not to execute with optimization and which variable is the "cost function" 
that is being minimized or maximized (the dependent variable). These variables are: 


TOLOPT 

NJOPT 

NVOPT 


Criteria of convergence on the dependent variable. The default is 0.0002. 

Component number which indicates the location of the dependent variable (if 
zero, the dependent variable is a performance property, not a DATOUT 
parameter). 

To turn off the optimization, NVOPT must be set to zero or all OPTVs must be 
turned off. If NVOPT is positive, the optimum is a minimum; if NVOPT is 
negative, the optimum is a maximum. 

If NJOPT equals zero: 

A value of 1 to 17 indicating which performance property is the dependent 
variable (see list of available performance properties in description of CNTL). 
If NJOPT not equal to zero: 

A value of 1 to 9 indicating which DATOUT of component number NJOPT is 
the dependent variable. 
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OUTPUT 

There is no DATOUT output array for OPTV components. 


Optimization usage notes: 

The optimal variable component (OPTV) allows the user to maximize or minimize a “cost 
function" (i.e. f the performance property or component output variable being optimized). An 
OPTV that is turned on (SPEC(9,N)=1) tells the program what variable (or variables if the user 
has multiple OPTV s that are turned on) the program can vary to minimize or maximize a “cost 
function . OPTVs are only used when optimization has been turned on (using variable 

NVOPT). Therefore optimization only takes place if NVOPT ± 0, and at least one OPTV is 
turned on. 

As an example of the use of an OPTV, let us assume that we have marched to Mach 1 4 
at 40000 feet and then throttled back to 50 percent F/Wa (See the previous example on 
marching in the description of CNTLs). We can now set SPEC(1,20)=1 to hold the FAVa at the 
present value. If we want to minimize the TSFC while holding F/Wa constant and varying the 
turbine inlet temperature, we would do the following: 

Assume that component 5 is the main burner. To perform the desired optimization a 
new component would be input as follows: 

KONFIG(l,21)=‘OPTV ,0,0,5,0,SPEC(1,21)=0,0,0,4,0,0,0,0,0, 

* n P u ^ record says that DATINP(4) (the burner outlet temperature) of component 5 
(the main burner) is the independent variable. There is no minimum or maximum value and 
since SPEC(9,21)=0, it is off. 

to™ N ° W to activate the optimizer we set SPEC(9,21)=1 and NVOPT=5 to minimize the 


Debugging: The variable DEBUG can be used to get intermediate output during 

optimization. If DEBUG=1, the program will print to the output dataset the following values: 


New value of the cost (unction 
' Original value of the coat function’ 

2 ) New value of the penalty function 
Original value of the penalty function’ 

( =1, if no variables have exceeded their user -defined limits), 

o) New value of beat value of the coat function 
Original value of the coat function ’ 

4) New vnhiea of the optimisation variable. 

Original valuea of the optimisation variables * 
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INPUTS 

The inputs at the design point are 
KONFIGU.N) 'Ll MV or 4HLIMV or 14 


KONFIGK2,N) 0 
KONFIG(3,N) 0 
KONFIG(4,N) 0 
KONFIGK5.N) 0 

SPECd^T) Blank 

SPEC(2,N) Minimum allowable value. 

SPEC(3,N) Maximum allowable value. 

SPEC(4,N) Type of limit variable being specified. 

= STAF , 4HSTAP , or 100 if the limit variable is a station property. 

= DOUT, 4HDOUT, or 200 if the limit variable is a DATOUT. 

= TERF , 4HPERF, or 400 if the limit variable is a performance property. 
SPEC(5,N) =Flow station containing the limit variable if it is a station property. 

^Component number containing the limit variable if it is a DATOUT. 

=0 if the limit variable is a performance property. 

SPEC(6,N) Component number, flow station number, or blank. 

SPEC(7,N) Blank 
SPEC(8,N) Blank 

SPEC(9,N) On-off switch, 1 = On, 0 = Off. 


OUTPUT 

There is no DATOUT output array for LIMV components. 

Variable limit usage notes: 

When a limit defined by the LIMV variable is exceeded, a warning message will be 
printed on the output file. If optimization is in effect, the criteria of merit, or the "cost function", 
will be penalized by a penalty function to drive the user away from exceeding the minimum and 
maximum limits. 
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input 

K0NFIG(1,N) 'SKED' or 4HSKED or 15 

KONFIG(2,N) Component number containing the dependent (scheduled) variable. 

KONFIGK3,N) Flow station number or component number containing the first independent 
variable on the map. X 

KONFIG(4,N) Flow station number or component number containing the second independent 
variable on the map. Y 

KONFICK5 ,N) Flow station number or component number containing the third independent 
variable on the map. Z 

SPEC(1,N) =0 turn schedule off, 

=1 turn schedule on. 

SPEC(2,N) SPEC number of the dependent variable of component KONFIGK2.N). 
SPEC(3,N) Scale factor on dependent variable (usually input as 1). 

SPEC(4,N) Table reference number containing the schedule map. 

SPEC(5,N) Location of first independent variable. 

I- 8 = STAP(1) - STAP(8) of flow station X (KONFIG(3,N)). 

II- 19 = DATOUT(l,N) DATOUT(9,N) of component X (KONFIGK3.N)). 

21-35 = SPECU.N) - SPEC(15,N) of component X (KONFIGK3,N)). 

SPEC(6,N) Scale factor on SPEC(5,N) (usually input as 1). 

SPEC(7,N) Location of second independent variable. 

I- 8 = STAP(l) - STAP(8) of flow station Y (KONFIGK4.N)). 

II- 19 = DATOUT(l,N) - DATOUT(9,N) of component Y (KONFIG(4,N)). 

21-35 = SPEC(l.N) - SPEC(15,N) of component Y (KONFIG(4,N)). 

SPEC(8,N) Scale factor on SPEC(7,N) (usually input as 1). 

SPEC(9,N) location of 3rd independent variable. 

I- 8 = STAP(1) - STAP(8) of flow station Z (KONFIG(5,N)). 

II- 1.9 = DATOUTU.N) - DATOUT(9,N) of component Z (KONFIG(5,N)). 

21-35 = SPEC(1,N) - SPEC(15,N) of component Z (KONFIG(5,N)). 

SPEC(10,N) Scale factor on SPEC(9,N) (usually input as 1). 

SPEC(11,N)- 

through- 

SPEC( 15, N) Blank 

OUTPUT 







Variable scheduling usage notes: 


to three variables (valued ** specifie<i by a ma P that is a function of up 

SPEC of the first 11 types *" UWr to dc8 ^ i 

three independent variables. TliehideDendelJt!^ components) to be a function of up to 
below. The value of the SPEC is determined bv internal* f ™ 8pecified b y u *«r as indicated 

format is the same as for other map tables.^ ** terpolatm * a U8er auppUed map. The map 


the tem^re^^^ 7 ** pro *™ to - 

^ 


— ui; 

the compressor pressure ratio (DATOUT 9) of compressor 3 . 


the mass flow rate (station property 1 see Tabl* v • L 

number 22 , and X ’ lab,e 3) en tenng burner 11 at flow station 


the rpm (SPEC 1) of shaft 31. 


map format, thrt ha, burner 

■Pm. The table reference nun.be, i, 3001 and the fl ° w - *"« 

ill 1 . 


follows: 


wui^vii Vni 09 . 

Ah «al. factor, are 1,. end the echedul. i, turned B. schedule would be Input as 


KONFIG( 1 ,39)='SKED’, 1 1J .22 31 SPPT/t io\_i , , „ 

.u*J,Z4,3i,SFEC(l,39)=l,4,1.0,900l, 19,1.0,3,1.0,21,1.0, 


Hi 


D. 
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INPUT 

K0NPIG(1^) 

K0NFIG(2,N) 

K0NFIGK3.N) 

K0NFIG(4,N) 

K0NFIG(5,N) 


'VCNTor 4HVCNT or 16 
0 
0 
0 


SPEC(1,N) 


SPEC(2,N) 


New value for the SPEC variable that conditional controls will reset if 
tnggered. ^/ conditional control is being used to turn a variable control on or 
off a value of RPECU N) less than or equal to zero will turn the variable 

onnt ' m ? U r ° u fSPEC(1 > N) * reator than ze ™ will turn the variable 
control on. Note: if the tolerance of a variable control is input as 1, the 

program will use the default of 0.001 (or 0.0001 if NVOPT £ 1) 

Component number for the variable that the conditional control will reset 


_ n . , ^ * thra * SPECs are used se t up a variable that will trigger the conditional 

,, ‘ . 3 value of this variable is monitored in relation to the trigger value to determine if 

the conditional control will reset the desired SPEC variable. This va^able is ref^rrln TTic 
tr^r variable, and is not to be confused with the SPEC variable that the conditional control 6 

SPEC(3,N) The type of trigger variable . 

•nnrni 100 iftrigger variabIe is & station property. 

'DP^^wniMp 7, ' trigger variab,e is a data output (DATOUT). 

PE op ISipp ° r 3 .°® lf . tngger variab le is a data input (SPEC). 

RPFPraxn o RF ’ 4HPERF > or 40° if trigger variable is n performance property. 

C ’ Station property number of trigger variable if SPEC(3,N)='STAP ' 

Data output number of trigger variable if SPEC(3,N)=’DOUT\ 

Data input SPEC number of trigger variable if SPEC(3,N)='DINP’. 
Performance property number of trigger variable if SPEC(3,N)=’PERF. 

,N) Flow station number of trigger variable if SPEC(3,N)='STAP\ 

Component number of trigger variable if SPEC(3,N)='DOUT. 

Component number of trigger variable if SPEC(3 N)='DINP’ 

0 if SPEC(3,N)='PERF. ’ 


SPEC(4,N) 


SPEC(6,N) 


SPEC(6,N) 


SPEC(7,N) 


Trigger value, the limit that has to be exceeded (either greater than or less 

Mm^TOnt^PECX2)! n ^ t *° na ^ COntr<>1 reS6t SPEC Variable SPEC(12) ° f 

This SPEC is only necessary if the user defines more than one conditional 
control to reset one particular variable. conditional 

rs°PEmVvi di r ° nal C ° nt ?L^' n0t re8et the desired SPEC variable 
«PEC 12, N) of component SPEC(2,N)) more than once per case. This prevents 

duri^ u?S2l? ^c n a a ]e COntr ° l8 fr ° m re>etting ° ne SPEC Variab,C back and forth 

** dnlml SPEC v ' ri,b " tvery «- «* 


0 



T 


SPEC(8,N) 

SPEC(9,N) 


SPEC(10,N) 

SPEC(11,N) 


SPEC(12,N) 

SPEC(13,N) 

SPEC(14,N) 


SPEC(15,N) 


=0 Do not reset SPEC(12,N) of component SPEC(2,N) at end of case. 

=1 Reset SPEC(12,N) of component SPEC(2,N) to the value it had at the 
beginning of the case before processing the next case. 

=0 Conditional control is nonactive. 

=1 Set SPEC(12,N) of component SPEC(2,N) to SPEC(1,N) if triggering 
variable is greater than the triggering value, SPEC(6,N). 

=-l Set SPEC(12,N) of component SPEC(2,N) to SPEC(1,N) if triggering 
variable is less than the triggering value, SPEC(6,N). 

Blank 

This input variable is only used if the conditional control is resetting 
SPEC(9,N) of a variable control. 

=0 Do not reset the independent value of the variable control. 

=1 Reset the independent value of the variable control to the value it had at 
the beginning of the case before processing the next case. 

Number of the SPEC for component SPEC(2,N) that the conditional control 
will reset. Default is 9. 

Blank 

Blank (If the conditional control is used to reset SPEC(9,N) of a variable 
control, this location is used to store the original value of the independent 
variable of the variable control.) 

Blank (used to store the original value of the SPEC variable that the 
conditional control resets, before any conditional controls are activated). 


OUTPUT 

There are no DATOUTs for VCNTs. 


Conditional control usage notes: 


Conditional controls (VCNT) allow the users to define variables that the program will 
monitor. If a variable exceeds a preset, user-defined limit, the conditional control will reset 
another user-defined variable. This user-defined variable could be one va riab le input (SPEC) 
for any component, such as: 1) resetting burner temperature or bypass ratios; 2) turning 
variable schedules on or off; or 3) turning variable controls on or off. 

The uses for conditional controls are best explained by examples. For the first use of 
conditional controls, suppose the user wants to turn on burner augmentation at Mach numbers 
above 1.1, but does not want to input the variables to do that directly. The user just sets up a 
conditional control to reset the augmentor temperature to the desired level (turn on 
augmentation) if the Mach number is greater than 1.1 and another conditional control to reset 
the augmentor temperature to zero (turn off augmentation), if the Mach number is less than 
1.1. For a second example, which involves the second and third uses of conditional controls, 
suppose that a user has an inlet airflow schedule to set engine airflow versus Mach number 
The user is throttling the engine thrust back at various Mach numbers, but can only throttle 
the engine so far at a constant airflow. Then the airflow must be allowed to be changed or the 
cycle will not converge. The user would set up conditional controls to watch one cycle 
pwameter, such as compressor surge margin, or burner temperature to turn off the airflow 
schedule (SPEC(1,N) of the correct variable schedule component) and maybe turn on a variable 
control to vary airflow to meet the desired engine performance level. 



78 


One conditional control can only be used to reset one SPEC variable to one new value, or 
to turn one variable schedule or variable control either on or off. For example, to turn a 
variable control on at one cycle condition, but turn it back off during the same case requires 2 
conditional controls. At the end of that case, conditional controls can be configured to reset all 
variables that they have changed. The independent variables of variable controls, activated by 
conditional controls can also be reset to their original values. The default option is to leave all 
SPEC variables and variable controls in their present states and values at the end of that case. 

Conditional controls are handled in the following manner. Before the program actually 
starts to iterate on the solution, it checks all active conditional controls whose trigger variable is 
a data input (SPEC(3 t N)= , DINF). if a data input has already exceeded its limit, the program 
will reset the specified SPEC variable and write to the output dataset that it has reset a SPEC 
variable. Then the program starts execution of the case. During the case, the program reaches 
a solution. The program will check each active conditional control to see if it has exceeded the 
specified limit for the specified trigger variable. If a limit has been exceeded for any conditional 
control, the program will reset the specified SPEC variable and write to the output dataset that 
it has reset a SPEC variable. The program will check all active conditional controls. Even if 
only one conditional control has reset a SPEC variable, the program will re-execute the case To 
prevent the program from entering an infinite loop, it will only re-execute a case a maximum of 
10 times with any variables, variable schedules, or variable controls changed by conditional 
controls. At the end of the case, if SPEC(8,N) of the conditional control is equal to 1 the 
conditional control will restore the original value of the SPEC variable. Otherwise the SPEC 
variable will stay at its current setting. If the conditional control was turning a variable control 
either on or off, setting SPEC(11,N) of the conditional control to 1 would cause the conditional 
control to reset the independent variable of the variable control to its value from the beginning 

of the case Otherwise, the independent variable of the variable control will also remain at its 
current value. 
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10.3 TABLE DATA TNPTTTS 


Consider the table data as three dimensional, composed of a series of planes with each 
Plane assigned a value called Z Then, on each Z plane, the dependent 

tablels oft^ refem*dtoM^ma|X <a * )ac ' Ma **•> «* Y **>" »• tari' 

Each table has the following input format: 

Record 1 Table Reference Number (Integer, columns 2-5). 

Table Identification Label, (Alpha-numeric, columns 6-75). 

Z - Identifier (Alpha-numeric, columns 1-4). 

NZ - Number of Z values (Integer, columns 6-7). 

Z - Variable values, 7F10., beginning in column 11. 

If nMded, extra card, follow 10X.7F10. format Z value. must be in .Mending order. 

Y - Identifier (Alpha-numeric, columns 1-4). 

NY - Number of Y values (Integer, columns 6-7). 

Y - Variable values, 7F10., beginning in column 11 

If needed, extra cards follow 10X.7F10. format. Y values must be in ascending 


Record 2 


Record 3 


Record 4 


X - Identifier (Alpha-numeric, columns 1-4). 

NX - Number of X values (Integer, columns 6-7). 

X - Variable values, 7F10., beginning in column 11 

If neod'd, extra card, follow 10X.7F10. format. X values must be in .Mending 


Record 5 


F(X,Y,Z) - Identifier (Alpha-numeric, columns 1-4). 

NX - Number of X values (Integer, columns 6-7). 

F(X,Y,Z) - Variable values, 7F10., beginning in column 11 
If needed, extra cards follow 10X.7F10. format. 

These values correspond to the values in the X identifier record. 

Last Record-The last record must consist of the 3 characters EOT in columns 1-3. 


Bfistrictinna 

1) NZ, NY, NX may not be blank or zero. 

2) Maximum of 70 tables. 

3) NZ, NY, NX are limited to 100 or less. 

4) a zero or blank table reference number will halt the read. 

5> * reS ““ tab ' < ’ ,tora « e “ co,d ‘ *• «”>* Subroutine 
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^mKinoT?* Ch /u lCal - dl f S0Ciat10 " option enables Program to handle just about any 

combination of chemical species for which thermodynamic data is available. (The default 

MeTI^ T ^ the NNEP89 pro *™“ is preset ™th properties for air ^nd JP4 
fuel.) 1116 dissociation option uses a rewritten version of the computer program for calculation 

properties. ° hemiCal E< * uillbnum Compositions (CEC) (ref. 16), to calculate^hermodynamic 
Morri ,. *T°. use th e chemical equilibrium option, the user must set ICEC=1 in the global &D 

Sw k ^ p cfc °r? r !*r m to -i* *• cec Jisst £> 

ioiiows. me input must contain information about the reactanta if 

CECN^m” ? n A may ,^. nt f' n J^ er ™ 0 ^y nem ' c data ®^^®sii^ e d^^ng^pt^iigi , n a'wparate 11 

ss**af« *? *■ ■<«*«— “ 

user , , An “ •*» 

release of ft, program. The er.ci format 

thermeH nary ^erroodynamic daUset ^readVe^is^the o^^np^reejuired to 

amounts of each reacted b^i^fo^eL 8 The FARRAY OA^^Y^n^ays ra 
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Order 

Contents 

Format 

Column 

First line 

REACTANTS 

3A4 

1-9 

Next lines 

One input line for each reactant species (maximum 15). Air 
must be the first reactant Each input line contains: 




(1) Atomic symbols and formula numbers (maximum 5 sets) a 

5(A2,F7.5) 

1-45 


(2) Blank 


46-53 


(3) Enthalpy or heat of combustion^ 

F9.5 

54-62 


(4) State: S, L, or G for solid, liquid, or gas, respectively 

A1 

63 


(5) Temperature, in Kelvin, associated with enthalpy in (3) 

F8.5 

64-70 


(6) Blank, J, or 

A1 

71 


(7) Blank 


72 


(8) Combustion efficiency (default is 1.0) 

F8.5 

73-80 

Last line 

Blank 




a. Program will calculate enthalpy or internal energy (3) for 
species in the thermodynamic data set at the temperature 
(5) if zeros are in columns 37-38. 




b. Enthalpy if column 71 is blank or J. Units are cal/g-mole if 
blank or J/formula wt. if J. Heat of combustion is in 
BTU/lb if column 71 is B. 




Table 8. Reactant Input 
10.4.1 Using The CEC Debugging Parameters 

The CEC program has parameters that can be used to debug the operation of the 
NNEP89 subroutines, the CEC subroutines, or both. These parameters are IDBUG, NDBUG, 
LDBUG, and IDEBUG. To use these parameters, they must be set to a non-zero value and 
input in the CEC &INPT2 Namelist. The values and information printed out to the output 
dataset for IDBUG are listed in Table 9. The values and output for the other parameters are 
listed in reference 4. IDBUG is the major debug parameter for the NNEP89 code. NDBUG and 
LDBUG are for debugging the THERME subroutine, the interface between the NNEP89 code 
and the CEC subroutines. IDEBUG is a CEC debug parameter and is specifically for debugging 
the CEC subroutines. Since the CEC subroutines are called many times during NNEP89 
execution, if IDEBUG=1, a lot of IDEBUG information will be produced (At least 100 times more 
than a "normal" output). 
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IDBUG 


=0 


No debug output 



=1 

Subroutine names, component and flow station numbers 


=2 

Output for IDBUG=1 

Mass flows for CEC subroutines calls 


=3 

Output for IDBUG=2 

Convergence information for burners and gas generators 


Table 9. Debug Output for Variable IDBUG 



The CE?LTf^ n £,T '"J'r* 8 and U the s,me “ Sample 4 in Section 12.4. 

nnHnn Th^ V ^ ? i , 4 " descnb « d h «™ >■> detail to help explain the use of the CEC 
option. The entire input file is given in Section 12.4. 


© 



0 

Figure 8. Block Diagram of Air-Turbo Ramjet 
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J) 








#nmn , The ^ U8 * 8 "standard air" (assumed to be 75.46% nitrogen, 23.19% oxygen and 1 35% 
argon by weight or N x .5 606 O i4198 Ar .0093 by atomic formula). The gas generator^ a 

nuxtore ofliquid and gaseous methane (CH 4 ) for the fuel and liquid and gaseous oxygen (Oo) 

The ® eC0 . n J bu ™er (augmented operates at stoichiometric conditions using a 
muttare o. unsymmetncal dimethylhydraxine (atomic formula C 2 H 8 N 2 ), a typical jet fuel (Sith 

a hydrogen to carbon ratio of 0.163, atomic formula CH L9423 , roughly JP4) and a second 

pineal jet fuel (with a hydrogen to carbon ratio of 0.161, atomic formula CH, 918 «, roughly 

JP5). The mass ratio ofliquid to gaseous methane in the gas generator fuel stream is 6 to 1 
o^dixer stream is 82% liquid oxygen and lf% 

frac^n of unsymmetncal dimethylhydraxine in the second burner is 0.40, forWiUs0.*0 and 

***• jh?r for ^* ***»“ JSSS anSr^^ouTd h» e 

< K«K5;gS , l f“ WOuld •• r “"“‘ for 

the user wants the cycle to use nif or no? C -n, C t “ le ?’ ^ the first mutant , whether 

efficiency end h»tii T* h ‘ S 8 reatta " t H>«ies. The conibustion 

Figur^ below Iheiewecohimn^Iiimhe 6 reartjnt mput for the esample case is shown in 

m :xXtT4^ z? zzes: e are ^ 

no more CEC input, so the octant input i, followed by a blanked » -e'nd" 1 ’ 
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1 2 3 
12345678901234567890123456789012345 

REACTANTS 

N 1.5606 O 0.4198 AR .0098 
C 1.0000 H 4. 

Cl. H 4. 

0 2 . 

O 2. 

c 2. H 8. N 2. 

Cl. H 1.9423 
Cl. H 1.9185 


€78901234567890123456789012345678901234567890 


00 

00 


100.0 

1 

bo 

00 

N> 

100.0 

100.0 

-21390. 

100.0 

100.0 

-3102 

100.0 

11900. 

100.0 

18300. 

100.0 

18600. 


G 298.15 O 
L 111.66 F 
G 298.15 F 
L 90.18 O 
G 298.15 0 
L 298.15 F 
1* 296.15BF 0.98 
1> 298.15BF 0.98 


END 


Figure 9. Example REACTANT Input Information 

each specie used and where iUrtM be^d* The^ARRAY^ **** P T ogra,a the rel *tive amounts of 
amount of each specie to be ad Ted Ts u USed f° tel1 the relativ « 

and gas generators. The OARRAY array is used to tell th* r f T ^ " ° n y reqmred for owners 
added as the oxidizer, and is only reauireH for «« th . C re a ^ ve amount each specie to be 
arrays are very flexible a ^ FARRAY and 0ARRAY 
for another. The same reactant sDeeie “ M . AW !^ Y f J °” e «° m Ponent can be in the FARRAY 
in thegas generator component, if necessary. FARRAY and O AKRA^ o'!? 0ARRAY &rray 
(ij). The j index corresponds with the component™™ ber of th^h^ 2dlTnension »* arrays, 

index contains pairs of reactant numhor nnH - l * number of the burner or gas generator. The i 

relative amount, second reactant number its refa^VT 88 '^ 11 “ i ? mt 88 8 reactant number, its 
reactants can be included inTach F^Eay or SSj r^’ “I" 1 80 ° n Up 10 6 different 
each reactant specie as it is read i nSlSS A is to 

methane is 2, gaseous methane is ^ i given example, air is 1, liquid 

relative mass indices The program will sum Inrf ° 6 n ^. mbe ™ or Actions can be used for the 
for the e^,e P XF^ n rSrh^r aliVe — “*"* 

FARRAY (1,5) =2, 6, 3,1, 

OARRAY (1, 5) =4, 82, 5, 18, 

FARRAY (1,6)- 6, .4, 7,. 3, 8,. 3, 


I1 M . 3 IMnp The rnmir. ru,.;- , Vi , hn n t rvr 

« withe* tSS£ Ferenc 

&D Namelist input. In the CEC 4INPTO N*tlv ^ \ done by 8etting ICEC ^ 1 in the global 
the desired value from Table 9 This ontinn h il In ? ut ’ vari ables ICEC=0 and IDBUG to 
error, in input b . M bwn f °“" d » >- very „„f u | f„ r d , tennining 

the input datawt would look like the following: P " dunn * executlon - The flrot few lines of 
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TITLE CARD 

CD ICEC-1 (and other global «D Namelist variables) 
(END 

NAMELIST 

(INPT2 IDBUG~1 , ICEOO, (END 
END 

(D (KONFIG and SPEC input) 


(END 


a 
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12.0 SAMPLE INPUT AND OUTPUT 


CO n^ OT ,r^ e i 0 < , 1,0V T g * 8eCti0nS ] l Bt th * input ffles for five Cerent engine configurations The 
corresponds to the third umnl. innill id alsA 




t! 


3 


J 


•i 

J 



12.1 SAMPLE 1 ■ TWO-SP OOL MTXRr>-FT i QW TtlRROFAN 

1 " \ 2 ; 8 rJ.’ mixed ' flow turbofan using the default thermodynamics routines. 
.f 89 "V draw a block dia Br am of the engine components (DRAW=T), and print information 
(number of iterations^altitude, Mach number, inlet recovery, engine airflow, gross thrust, fuel 
tlow, net thrust, and TSFC) to the terminal as the case is executing (1TERM=2) A block 

diagram of this engine in shown in Figure 10. This case also includes input that is commented 
out using r and */. 

SAMPLE MIXED FLOW TURBOFAN 
*D I TERM- 2, DRAW-T , LONG-T, TABLES-T, PINPUT-T fcEND 
*D MODE-1 

KOHFIG' , 1)-4HINLT, 1,0, 2,0, SPEC (1,1) -100. 0, O . , 0 . , 0 . , 0 . , 

0.,0.,0.,0.,0.,0.,18.00, ' ' 

fn^ IG 1 <1 A 2 !^5 COMP ' 2 ' 0 ' 3 ' 0 ' SPECa ' 2 >" 1 - 500 ' 0 - .1-0. 1001., l.o, 

1002.. 1.0. 1003. .1.0. 0..0. .0.8500.3.0.1.0, 

KONFIG (1 , 3) -4HSPLT, 3 , 0, 4 , 9, SPEC (1 , 3) -1 . 0, 0 . 2000E-01 0 2000E-01 

SPEC<1 ' 5 > -O-S&OOE-Ol . °. 30° 0 ,0. . 3000. , 0 . 9900. 

fAnr lG , <1 A 6) “ 4HraRB ' 6 ' 12 ' 7 '°' SPEC < 1 -6)-3.500, 0.7500, 1.0 1007 1 0 

1008 ., 1 . 0, 1 . 0, 0 . 8000, 1 . 0, 0. 9000, 5000 ., 1 0 *' ' ' 

^mft IG i 1 n 7 i“n H T U A B ( 7 A 12 ' 8 '°' SPEC(1 ' 7 > , * 2 ' 5 0 0 , 0 -2500, 1.0, 1009., 1.0 

1010. ,1.0, 1.0, 1.0, 1.0,0. 9000, 5000. ,1.0 ' ' 

vn2^^M!- 4HMIXR ' 8 ' 9 ' 10 '°' SPEC <l'8>-0-.0., 0.4000, 1.0, 

? 9 ”n HN0ZZ ' 10 ' 0 ' 11 ' °' SPEC(1 ' d> "°-' 0 - 9800 '°-.0., 0.9750, 

1 . u , u . 1,0, 

KONFIG(l, 10) -4HLOAD, 0, 0, 0 , 0 , SPEC ( 1 , 10) — 200 . 0 , 

irnwrir !i ’ i J ! -4HSHFT, 4, 6, 10, 0, SPEC ( 1 , 11 ) -8000 . , 8*1, 

1.12 -4 HSH FT, 2 , 7 ,°, 0 ,SP ec ( 1 , 1 2) - 6 ° 0 °., 8*1 ,1.0,0., 0., 

KONFIG (1 , 13) -4HCNTL, SPCNTL (1, 13)-1. 0,7.0, 4HSTAP 8 0 0 0 1 

KONFIG(l,14)-4HCNTL,SPCNTL(i:i4Ul.o:6:o:4HiTAP 8 0 7 6 1 ' 

KONFIG (1, 15) -4HCNTL, SPCNTL (1, IS) *1.0, 4.0, 4HSTAP* 8 " 0 # 6 # 0* 1* 1 1 1 7S 
KONFIG (1, 16) -4HCNTL. SPCNTL 1 16 - 1.0 3 0 4HSTAP 8 0 4 o'l' ' 

KONFIG (1, 17) -4HCNTL, SPCNTI. 1 17 -1.0 1 0 4HSTAP 8 0 2 o'l' 

KONFIG (1 , 18) -4HCNTL, SPCNTL (l)l8)—l. 0^2.0 4HDOUT 808011 1 5 ^ 

KONFIG(l,19).4HCNTL,SPCNTL(i:i9!-1.0:n:00,5SS5{j?.i:^?i 0 i ' 

KONFIGil, 2 0,-4HCNTL. SPCNTL (1, 20) -1 . 0, 1 2 . 00 , 4HDOUT, 8 . 0^ 12," 0," l[ 

/* ^ U °- 3 °'^ TP “ 5000., ETAR-1. 0000, MODE- 1, SEND 

0 • 30, ALTP- 10000., ETAR-1. 0000, MODE- 1 SEND 
4D MACH- 0.30, ALTP- 15000 ., ETAR-1 . 0000, MODE- 1, SEND 
0.30, ALTP- 20000., ETAR-1. 0000, MODE- 1, SEND 
0.30, ALTP- 25000., ETAR-1. 0000, MODE- 1, SEND */ 

0.50, ALTP- 15000. , ETAR-1 .0000, MODE- 1, SEND 
0.70, ALTP- 25000., ETAR-1. 0000, MODE- 1, SEND 
0.80,ALT p - 40000., ETAR-1. 0000, MODE- 1, SEND 
1.20, ALTP- 40000. , ETAR-0. 9915, MODE- 1, SEND 
1.40,ALT p - 40000., ETAR-0. 9782, MODE- 1, SEND 
1.80, ALTP- 40000., ETAR-0. 9445, MODE- 1, SEND 


SD MACH- 
SD MACH- 
40 MACH- 
40 MACH- 
40 MACH- 
40 MACH- 
40 MACH- 
40 MACH- 
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12.2 SAMPLE 2 - SEPARATE FLOW TURBOFAN USINO OFO 


Sample 2 is a 1-spool, separate-flow turbofan (with regenerator) using the CEC 
thermodynamics routines (ICEC=1) and automatic controls (ACTL=1) to take care of all flow 
and work errors. The bypass air is being heated by the core flow before the core flow exits 
through a nozzle. The fuel is a "standard" JP fuel with a heating value of 18300.BTU/lb and a 
burner efficiency of 0.99. A block diagram of this engine in shown in Figure 11. 


TEST CASE: HTEX SANFORD NEW 

4D ACTL-1, ICEC-1, LONG-F, NCODE-1, 
REACTANTS 

N 1.5606 0 0.4198 AR .0098 
C 0.5245 H 1.0000 


DRAW-T, NMODES-1, SEND 

100. 0.0 G 298.15 O 

100. -4628. 066L 298.15 F 


NAMELIST 

4INPT2 ICEC-1, 4 END 
END 

4D MODE-1, 


KONFIG (1 , 1 ) -4HINLT, 1 
KONFIG (1 , 2 ) — 4HCOMP, 2 
1,0,0, .85, 12.0,1, 


,0,2,0, SPEC (1 , 1) -5*0, .97,5*0,14.4,0,100, 

,0,3, 12, SPEC (1, 2) -1 . 3, . 001, 1, 1004, 1, 1005, 1, 1006, 


KONFIG (1, 3) -4HSPLT, 3,0 
KONFIG (1,4) -4HDUCT, 4 , 0 
FARRAY (1,4)— 2., 1., 


,9, 4, SPEC (1,3) -20. ,0,0, 

,5,0, SPEC (1,4) -.05, 0,0, 3200. 


.99, 18300, 


1? N 91?500o!’! HTURB ' 5 ' 12 ' 6 ' 0 ' SPEC(1 ' 5) " 3 - 3 ' 1 ' 1 ' 1007 ' 1 ' 1008 ' 1 ' 1 '- 1 ' 


KONFIG(l, 6)-4HHTEX, 9, 6, 10, 7, SPEC(1, 6) -.05, . 03, 800 8 l 
KONFIG (1 , 7) — 4HNOZZ, 7 , 0, 8, 0, SPEC (1,7)— 0,1,6 jo,0?9850 11 0 1 
KONFIG(l,8)-4HSHFT,2,5,0,0,SPEC(i;8)-9n.C ' ' ' 

4 END 16 (1 ’ ^ “4HNOZZ , 10,0,11,0, SPEC (1 , 9)— 0, 1 , 0, 0, °, 9850, 1 , 1 , °, 1, 


4D MACH- 
40 MACH- 


0.30, ALTP- 5000. ,ETAR-1. 0000, MODE- 1, SEND 
0.50, ALTP- 15000. ,ETAR-1. 0000, MODE- 1, SEND 


© 



Q 


Figure 11. Block Diagram of One-Spool, Separate-Flow Turbofan (with Regenerator) 
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12.3 SAMPLE 3 - ONE-SPOO L MIXED-FLOW TURBOFAN 


Sample 3. is a 1-spool, mixed-flow turbofan using the default thermodynamics routines 
and output of the convergence history (LONG=T). This case also uses the XNUM array to name 
the flow stations on the output, but XNUM does not work on systems that cannot handle 
Namelist Hollerith or quoted string input. Removing the XNUM array would not cause any 
problems executing this input case. This case also has the headings option on (DOUTHD=T). 
(Prints out headers on the output file to aid the user interpreting the output, highly 
recommended for novice users.) The program will draw a block diagram of the engine 
(DRAW— T) and will make full passes through the engine while it determines the error matrix 
(NCODE=-l). A block diagram of this engine in shown in Figure 12. 


TEST CASE: MIXED FLOW TURBOFAN SANFORD NEW 

‘D LONG-T, NCODE — 1, DOUTHD-T, DRAW-T, 4 END 
SO MODE-1, 

XNOM-4HINLT, 4HLPCI, 4HSPLI, 4HHPCI , 4HM3RN, 4HHPTI, 4HMXRI, 4HABRN, 
4HNOZI, 4HNOZO, 4HSECI, 4HSECO, 

KONFIGU, 1)-4HINLT, 1,0, 2,0, SPEC (1,1) -100. 0, 0 . , 0 . , 0 . , 0 . , 0 . , 0 . , 
0.,0. , 0 . , 0 . , 18.00 r 

?SS3^?, < S:?i:^ c S?f2 2 9?6 3 9S: SPEC<1 ' 2, " 1 - 60 '°-' :L ' 1001 ' 1 ' 1002 ' 1 ' 

KONFIG (1, 3) -4HSPLT, 3,0,4,12, SPEC(1, 3)— 1 . 75, 0. 02, 0. 02, 

FARRAY(l'5)l2 H l UCT ' 5,0 ' 6 ' 0 ' SPEC<1 ' 5> “ ,05 '°' 0 ' 3000 -'- 99 ' 18300 -' 

KONFIG <1, 6) -4HTURB, 6, 0, 7, 0, SPEC (1, 6) -3. 6, 1 . 0, 1 , 1007, 1 , 1008, 

1,1, 1,1, .90,5000,1, 

KONFIG (1, 7) -4HMIXR, 7,13,8,0, SPEC (1 , 7)— 0.,0., 0.40,1. ,0 0 

FAR^Y!i;8)r2 H l UCT ' 8 ' 0 ' 9 ' 0 ' SPEC<l ' 8> ‘' 05 ' 0 ' 0 ' 3600 -'- 99 ' 1 ® 300 -' 

KONFIG (1, 9) -4HNOZZ, 9,0, 10,0, SPEC (1 , 9) -0 . , 0 . 98, 0 . , 0 . , .98,1.,1., 

FARRAY(i;i0)l2 H l OCT ' 12 ' 0 ' 13 ' 0 ' SPEC<1 ' 10> “* 05 ' 0 ' 0 ' 0 ’' 0 - 99 ' 18300 -' 
KONFIG (1 , 11) -4HSHFT, 2,4, 6, 0, SPEC ( 1 , 1 1 ) -9*1 . 0, 

If) -4HCNTL, SPCNTL <1, 12 > -1, 1, 4HSTAP, 8, 2, 0, 0, 

KONFIG (1, 13) -4HCNTL,SPCNTL(1, 13) -1,4, 4HSTAP, 8,4,0 0 
KnNPTr!! ' !!! SPCNTL (1>14)— i'll, 4HDOUT ,8,11*0^0, 

! ’1 HCNTL ' SPCHTLd. 15) -1, 6, 4HSTAP, 8, 6, 0, 0, 

KONFIG (1 , 16) -4HCNTL, SPCNTL (1 , 1 6) -1 , 3, 4HDOUT, 8, 7 , 0, 0 
17 > -4 hc NTL, SPCNTL (1 , 1.7) -1 , 2 ; 4HDOU T ; 5* 2 ) 50 ,1 , 

(END 

! 8 ^J 9 'J? ) " 1 ' SPEC < 9 ' 13 )- 1 -SPEC(9,14)-l,SPEC(9,15)-l, 

SPEC (5,1) -.26, SPEC (9,1) -50, SPEC (6, l)-0, 

SPEC (9, 16)— 1, SPEC (9,17)— 1, (END 








12 A SAMPLE 4 - AIR-TURBO RAMJET TIRING rvr 


Sample 4. is a 1-spool, gas-generator Air-Turbo Ranriet, using the CEC thermodynamic 
routines (ICEC=1). NNEP89 will draw a block diagram of the iycle ©RAW=52Se 
information about the convergence history (LONG=T). The code will make a full pass through 
the cycle when it calculates the partial derivative matrix (NCODE=-l). The gas-generator 
operates on liquid and gaseous methane for the fuel and liquid and gaseous oxygen for oxidizer 
combustion occurs at equivalence ratios greater than stoichiometric. The cycle has a 
stoichiometric augmentcr fuelled with a mixture of unsymmetrical dimethylhydrazine (CjjHgNj) 

and 2 hydr^rW species similar to JP4 and J?5. The CEC subroutines will determine the 
enthalpies for the hydrocarbon species from the lower heating value input in that reactant line 
The enthalpy for gaseous methane and oxygen will be calculated by the program because of the 
00 in columns 37-38. This input file also contains Namelist input for the CEC program with 
the variable IDBUG. Since IDBUG=0, debugging is turned off, but setting IDBUG to 1, 2 or 3 
would print our additional information that could be very useful for finding errors when an 
input case will not execute. Caution is suggested when using the CEC debug parameters it 

increases program output by a factor of about 10 or more. A block diagram of this engine in 
shown i Figure 13. 


,n T j™ ? EM ST E t^ " SE ° F CEC ' GAS GENERATORS AND F ARRAYS 
4D ICEC-1, DRAW-T, LONG-T, NCODE— 1, (END 
REACTANTS 

’606 O C.4198 AR .0098 
H 4. 

H 4 . 00 


l.: 

l. 

1. 

2 . 
2 . 
2 . 
1 . 
1 . 


8. N 

1.9423 

1.9185 


2.0 


00 


1. 

100.0 

100.0 

100 . 

100 . 

100.0 

100.0 

100.0 


-28.2 

-21390. 


-3102. 

11900. 

18300. 

18600. 


NAMELIST 

(INPT2 IDBUG-0, (END 
END 

(D MODE-1, 

KONFIG (1,1) -1,1, 0,2,0, SPEC (1,1) -100, 4*0, 
KONFIG (1,2) -4, 2, 0,3,0, SPEC (1,2) -1.6, 0,1, 

1.0. 0. .85.2.0.1., 

KONFIG (1,3) -7, 3, 0,4, 5, SPEC (1, 3) -0.,0., 
KONFIG (1,4) -8, 4, 10, 7,0, S.?EC (1 , 4) -0 . , 0 . , 0 
*.ONFIG (1.5) -3, 5, 0,6,0, SPEC (1,5) -2960, .67 

FARRAY (1,5) -2, 6,3, 1, 

OARRAY (1, 5) -4, 82, 5, 18, 

KONFIG (1, 6) -2, 7, 0, 8, 0, SPEC (1, 6) -.00001. 0 
SPEC (11 , 6) -1 , 

FARRAY (1,6) -6, .4,7, .3,8, .3, 

KONFIGG, 7) -9, 8, 0,9,0, SPEC (1,7) -0,1, 0,0, 

KONFIG (1 , 8) -5, 6, 0,10,0, SPEC (1, 8) -3 . 6, 6 

1008.1.1.0. 1. .80.5000.1, 

KONFIG (1,9) -11, 2, 8, 0,0, SPEC (1,9) -9*1, 
KONFIG (1,10) -12, SPCNTL (1,10) -1.0, 1.0, 100 

KONFIG (1,11) -12, SPCNTL(1, 11 >-4.0,5.0,200 
KONFIG (1 , 12) -12 , SPCNTL (1 , 12) -3. 0,5.0, 100 

KONFIGG, 13) -12, SPCNTL (1,13) -1.0, 8. 0,200 
KONFIG (1,14) -12, SPCNTL (1,14) -1.0, 9. 0,200 
KONFIG G, 15) -i2, SPCNTL (i, i5) -4. 0,5.0, 100 

tD SPEC (9, JO) -1, SPEC (9, 11) -1 , SPEC <9, 12) - 
SPEC (5,1) -.26, SPEC (9,1) -50, SPEC (6, l)-0, 
(D SPEC (5, 1) -. 6, SPEC (9, 1) -10000, (END 


.15 

.66 


298 
111 

298.15 
90.18 
298.15 
298.15 F 
298.15BF 
298. 15BF 


.97,2*0, .1,2*0,14. 
1001,1,1002,1,1003, 


25, .95, 

184, 600. ,1.644, 


0. 0,0. '.95, 21500., 


•985,1,1,0,1, 
1 , 1007 . , 1 , 


, 8 . 0 , 2 . 0 , 0 ., . 0 , 0 . , 0 ., 

, 8 . 0 , 4 . 0 , 0 , 0 , 0 . , 0 . , 

, 8 . 0 . 6 . 0 , 0 ., . 0 , 0 . , 0 . , 
,8.0, 9.0,0., . 0, 0. , 0 . , 
,5.0,2.0,20., .0, 0. , 0. , 
,2.0,10.0,28.510,0,0,0, 

1,SPEC(9,13)-1, 

SEND 


0.9e 

0.98 


99 



tD SPEC (5, 1) -. 9, SPEC (9,1) -25000, tEND 
tD SPEC (5, 1) -1 . 2, SPEC (9, 1 ) **32200, tEND 
iD SPEC (5,1) -1.5, SPEC (9,1) -39500, tEND 
tD SPEC (5,1) -1.8, SPEC (9,1) -46800, tEND 
tD SPEC (5 , 1) -2 . 0, SPEC (9, 1) -50900, tEND 
tD SPEC (5, 1) -2 . 5, SPEC (9, 1) -55000, tEND 
tD SPEC (5, 1) -3 . 0, SPEC (9,1) -58400, tEND 





12.5 SAMPLE 5 - TANDEM F AN WITH MtlLTIMODER 


Sample 5. is a 2-spool, tandem fan engine, with multimode capabilities. Mode 1 is in 
series, mixed flow, mode 2 in parallel-flow, high bypass. The input is set to draw a block 
diagram of each mode (DRAW=T), and the mode to design the components is mode 1 
(MODESN=l). This is a 2 mode engine (NMODES=2) and will make full passes through the 
engine while it determines the error matrix (NCODE=-l). The data output headers option is on 
(DOUTHD=T). (Prints out headers on the output file to aid the user interpreting the output, 
highly recommended for novice users.) This case will print information (number of iterations, 
altitude, Mach number, inlet recovery, engine airflow, gross thrust, fuel flow, net thrust, and 
TSFC) to the terminal as the case is executing (ITERM=2). This engine uses XNUM to label the 
station property outputs, the XNUM array could be removed without changing the way the 
engine is operated. Since this engine may require more than 50 iterations to converge, the 
maximum number of iterations has been increased to 75 (MAXNIT=75). A block diagrams of 
this engine in modes 1 and 2 are shown in Figures 14 and 15, respectively. 


2 MODE TANDEM-FAN ENGINE 1-SERIES, MXDFLO 2-PARALLEL HY-BY 
SD NCODE— 1 , DRAW-T , LONG-F, NMODES-2 , MODESN-1 , DOUTHD-T, SEND 
SD MODE-1, ITERM-2. MAXNIT-75, 

XNUM-4HINTI , 4HDSPI , 4HFFNI , 4HFFNO, 4H2FNI , 4H2FNO, 4HHPCI , 4HHPCO, 4HHPTI 
4HLPTI, 4HMXMI , 4HMIXO, 4HNOZI , 4HNOZO, 4HDDTI , 4HDNZI , 4HDNZO, 4HCOOL, 4HBDUC, 
4HMXSI , 

KONFIG (1, 1)-4HINLT, 1,0, 2,0, SPEC (1,1) -300. 0, 0. , 0. , 0. , 0. , 1.0, 

KONFIG (1, 2) -4HCOMP, 3, 0,4,0, SPEC (1,2) -1.5,0., 1.0, 1001., 1.0, 

1002. , 1.0, 1003. , 1. 0, 00. 00, 0., 0.8800, 2. 800, 1.0, 

KONFIG(l, 3) -4HCOMP, 5, 0,6,0, SPEC (1,3) -1.5,0., 1.0, 1001., 1.0, 

1002. . 1.0. 1003. . 1. 0. 00. 00. 0.. 0.8800. 1.500. 0.95, 

KONFIG(1,4)-4HSPLT,6,0, 7, 19, SPEC (1,4 ) -0.7000, 0 . 2000D-01 , 0 . 2000D-01 , 

KONFIG (1, 5) -4HCOMP ,7,0, 8, 1 8, SPEC (1 , 5 ) -1 . 3 , 0 . 1 00QD00, 1 . 0, 1004 . , 1 . 0, 

1005. , 1 .0, 1006. , 1. 0, 0. , 0. , 0.8800, 7.200, 0. 935, 

KONFIG (1, 6) -4HDUCT, 8,0,9. 0, SPEC (1, 6) -0 . 6000D-01 , 0 . , 0 . , 3260. , 0. 9850 
0. 1830D05, 0. , 0. , C. , 0 . 1000D00, 

KONFIG (1,7) -4HSPLT, 2,0,3, 15, SPEC (1 , 7) -0 . 10000-01 , 

KONFIG (1, 8) -4HTURB, 9, 1 8 , 1 0 , 0 , SPEC (1 , 8) -3 . 500, 0 . 7000, 1 . 0, 1007, 1 . 0, 

1008, 1.0, 1.0, 0.5000, 1.0, 0.9000, 5000. ,1.0, 

KONFIG (1 , 9) -4HTURB, 10, 18, 11,0, SPEC (1,9) -2. 500, 0.3000, 1.0, 1009, 1.0, 

1010.1.0. 1 .0.0.5000. 1 .0. 0. 8800. 5000. . 1 . 0, 

KONFIG ( 1,11) -4HDUCT ,19,0,20,0,SPEC(1,11)-0.5000D-01,0. , 0 ., 400 . 0 , 0 . 9700 
0. 1830D05, 0. , 0. , 0. , 0.5000D-01, 

KONFIG(1,12)-4HMIXR,11,20,12,0,SPEC(1,12)-0.,0., 0.3000, 0.9500, 1.0,1 0 
KONFIG (1, 14) -4HDUCT , 12, 0, 13, 0, SPEC ( 1 , 14 ) -0 . 5000D-01 , 0 . , 0 . , 400 . 0, 0 . 9700 
0.1830005,0. , 0. , 0. , 0.5000D-01, 

KONFIG (1, 15) -4HNOZZ, 13, 0,14,0, SPEC (1, 15) -0 . , 1 . 0, 0 . , 0 . , 0 . 9800, 1.0,0.,0. 1 0 
KONFIG ( 1 , 16) -4HLOAD, 0,0, 0,0, SPEC (1,16)— 100.0, ' ' 

KONFIG(l, 17) -4HSHFT, 5, 8, 0,0, SPEC (1,17) -0.1400005,1. 0,1. 0,1. 0,1.0, 

1*0, 1.0, 1.0, 1.0, 

KONFIG (1, 18) -4HSHFT, 2,3, 16, 9, SPEC (1 , 1 8) -8000 . , 1 . 0, 1 . 0, 1 . 0, 1 . 0, 

1.0. 1.0. 1.0. 1.0, 

KONFIG (1,19) -4HDUCT ,15,0,16,0,SPEC(1,19)-0.5000D-01,0.,0.,400 0 0 9700 
0.1830D05, ' 

KONFIG (1, 21) -4HCNTL, SPCNTL (1,21) -1.0, 9.0, 4HSTAP, 8.0, 

13. 00. 0. .1.0. 1.0.0., 

KONFIG (1, 22) -4HCNTL, SPCNTL (1,22) -1 .0,8.0, 4HSTAP,8.0, 

10 . 00 . 0 .. 1 . 0 . 1 . 0 . 0 . , 

KONFIG (1, 23) -4HCNTL, SPCNTL (1,23) -1.0, 17 . 00, 4HDOUT, 8 . 0, 

17 .00. 0. . 1 .0. 0. . 0. , 

KONFIG (1 , 24 ) -4HCNTL, SPCNTL (1 , 24 ) -1 . 0, 1 8 . 00 , 4HDOUT, 8.0, 

18 . 00. 0 . . 1 . 0. 0 . . 0 . , 

KONFIG (1, 25) -4HCNTL, SPCNTL <1, 25) -1 .0, 2 .0, 4HDOUT, 5.0, 

2. 0. 25. 00.0. .1.050. 4.0, 

KONFIG (1, 26) -4HCNTL, SPCNTL (1, 26) -1 . 0, 3 . 0, 4HDOUT, 5.0, 

3. 0. 25. 00.0. .1.050. 4.0, 

KONFIG (1, 27) -4HCNTL, SPCNTL (1 , 27) -1 . 0, 5 . 0, 4HSTAP, 8.0, 
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9. 0. 0. .1.0. 1.050. 4.0, 

KONFXG (1 , 28) -4HCNTL, SPCNTL (1,28) >1.0, 4.0, 4HDOUT, 8.0, 

12.00. 0.. 1.0. 0..0., 

KONFIG (1,29) -4HCNTL , SPCNTL (1,29) -1.0, 1.0, 4HSTAP ,8.0, 

3.0. 0.. 1.0.0. .0., 

KONFIG (1,31) -4HCNTL, SPCNTL <1,31)-10.00,2.0, 4HSTAP ,8.0, 

5. 0. 0. .1.0.0. .40. 00, 

KONFIG (1, 32) -4HCNTL, SPCNTL (1 , 32) -10 . 00, 3 . 0, 4HSTAP, 8.0, 

7. 0. 0. .1.0.0. .40. 00, 

KONFIG (1, 37) -4HLIMV, 0, 0, 0, 0, SPEC (1,37) -0,0, 8800., 4HDOUT, 2., 18., 0,0,1, 

KONFIG (1, 40) -4HLIMV, 0, 0, 0, 0, SPEC (1,40)-0,0,1.100, 4HDOOT, 6. ,2. ,0,0,1, 

KONFIG (1, 45) -4HLIMV, 0, 0, 0, 0, SPEC (1,45) -0,0,1. 100, 4HDOUT, 6., 3., 0,0,1, 
KONFIG (1, 48) -4HLIMV, 0, 0 , 0, 0 , SPEC (1 , 48) -0, 1 . 050, 4 . 0, 4HDOUT, 4 . , 5. , 0, 0, 1, 
KONFIG (1, 49) -4HLIMV , 0, 0, 0, 0, SPEC (1, 49) -0, 0, 1 . 100, 4HDOUT, 6. ,5. ,0,0,1, 

KONFIG (1, 50) -4HLIMV, 0, 0, 0, 0, SPEC (1,50) -0,2. 0,80. 00, 4HDOUT, 5., 5., 0,0,1, 
KONFIG (1, 54) -4HOPTV, 0,0, 15,0, SPEC (1, 54) -0.,0.,0., 1.0,0., 

0.,0.,0.,1.0, 

KONFIG (1, 58) -4HCNTL, SPCNTL (1,58) -4.0, 6.0,4HDOUT, 6.0, 

3. 0. 1 . 0. 1 . 0. 0. . 0 . , 

KONFIG(l, 59) -4HDUCT, 4, 0,5,0, SPEC (1,59) -0., 

KONFIG(l, 60) -4HNOZZ, 16, 0,17,0, SPEC (1, 60) -0 . , 0 . lOOODOO, 0 . , 0 . , 0 . lOOODOO, 

1 .0. 1.0. 0. . 1 .0, 

4 END 

&D MODE-2, 

XNUM-4HINTI , 4HSPLI , 4HFFNI , 4HVLVI , 4H2FNI , 4H2FNO, 4HHPCI , 4HHPCO, 4HHPTI 
4HBNZl' 4HBNZO' 4H ' 4HN0ZI ' 4HN0Z0 ' 4HFFN0 ' 4HDNZI, 4HDNZO, 4HCOOL, 4HBDUC, 

KONFIG (1 , 1 ) -4HINLT, 1 , 0 , 2 , 0 , 

KONFIG (1, 2) -4HCOMP, 3, 0, 15, 0, 

KONFIG (1,3) -4HCOMP, 5,0, 6, 0, 

KONFIG (1 , 4) -4HSPLT, 6,0,7,19, 

KONFIG (1 , 5) -4HCOMP, 7,0,8, 18, 

KONFIG (1 , 6) -4HDUCT, 8,0, 9, 0, 

KONFIG (1, 7)-4HSPLT,2, 0, 4, 3, 

KONFIG (1, 8) -4HTURB, 9,18,10, 0, 

KONFIG (1 , 9) -4HTURB, 10, 18,11,0, 

KONFIG (1, 11) -4HDUCT, 19,0,20,0, 

KONFIG (1, 13) -4HNOZZ, 20, 0,21,0, SPEC (1, 13) -O . , 1 . 0, O . , 0 . , 0.9500, 

1 . 0 , 0 . , 0 . , 1 . 0 , 

KONFIG (1, 14) -4HDUCT, 11, 0, 13,0, 

KONFIG (1, 15) -4HNOZZ, 13, 0, 14, 0, 

KONFIG (1, 16) -4HLOAD, 0, 0,0, 0, 

KONFIG (1 , 17 ) -4HSHFT, 5 , 8 , 0 , 0 , 

KONFIG (1 , 18) -4HSHFT, 2,3,16, 9, 

KONFIG (1, 19) -4HDUCT, 14,0,16,0, 

! CO n F ^ G<1 ' 20) *' 4HNOZZ ' 16 ' 0 ' 17 ' 0 ' SPEC(1 ' 20) “ 0 -' 1 - 0 ' 0 -' 0 *' 0 - 9800 ' 

0. , 0. , 1 .0, 

KONFIG (1,22) -4HCNTL , SPCNTL ( 1 , 22 ) -1 . 0 , 8 . 0 , 4 HSTAP ,8.0, 

10. 00,0. ,1.0, 1.0,0., 

KONFIG (1, 23) -4HCNTL, SPCNTL (1,23) -1.0, 17. 00, 4HDOUT, 8.0, 

17.03.0. .1.0.0..0., 

KONFIG (1, 25) -4HCNTL, SPCNTL (1,25) -1.0, 2. 0,4HDOUT, 5.0, 

2. 0. 25. 00.0. .1.050. 4.0, 

KONFIG (1 ,26) -4HCNTL, SPCNTL (1, 26) -1.0, 3.0, 4HDOUT, 5 . 0, 

3. 0. 25. 00.0. .1.050. 4.0, 

KONFIG (1, 27 )-4HCNTL, SPCNTL (1,27) -1.0, 5. 0,4HSTAP, 8.0 

9. 0. 0. .1.0. 1.050. 4.0, 

KONFIG (1 , 29) -4HCNTL, SPCNTL (1 , 29) -1 . 0, 1 . 0, 4HSTAP, 8.0, 

3. 0, 0. , 1 . 0, 0 . , 0 . , 

KONFIG (1, 33) -4HCNTL, SPCNTL (1 , 33) -1 . 0, 18. 00, 4 HSTAP, 8.0, 

13 . 00, 0. , 0 . , 1 . 0, 0 . , 

KONFIG (1, 34) -4HCNTL, SPCNTL (1, 34) -1 . 0, 9.0, 4HDOUT, 8.0, 

18.00. 0.. 1.0. 0..0., 

KONFIG(l, 35) -4HCNTL, SPCNTL (1,35) -1.0, 4. 0,4HSTAP, 8.0, 

7.0, 0. , 1 .0, 0. , 0. , 

' oa! “ 4h ^ imv » °» °« °> °» S PEC (1,37) -0,0, 8800., 4HDOUT, 2. ,18. ,0,0,1., 
KONFIG (1, 39) -4HCNTL, SPCNTL (1, 39) -10 . 00, 3 . 0, 4HSTAP, 8 . 0, 

20.00. 0..0..0..40.00, 

KONFIG (1, 40) -4HLIKV, 0, 0,0,0, SPEC (1,40) -0,0,1. 100, 4HDOUT, 6., 2., 0.0 1 
KONFIG (1, 45) -4HLIMV, 0, 0, 0, 0. SPEC <1, 43) -0, 0 1 .100 4HDOOT 6 3 0 0 1 

KONFIG (1 , 48) -4HLIMV, 0, 0, 0, 0, SPEC (1 , 48) -0, 1 ! 050, 4. 0, 4HDOUT, 4?, 5?) o| 0," 1 . 
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KONFIG (1 , 4 9) -4HLIMV, 0, 0, 0, 0, SPEC (1,4 9) “0,0,1.100, 4HD0UT ,6. ,5. ,0,0,1., 
KONFIG (1 , 50) -4HLIMV, 0 , 0 , 0, 0, SPEC (1,50) “0,2.0,80.0, 4HDOUT ,5. ,5. ,0,0,1., 
KONFIG (1 , 51 ) -4HCNTL, SPCNTL (1 , 51) -10 . 00, 2 . 0, 4HSTAP ,8.0, 

16.00. 0. .0..0.. 40.00, 

KONFIG (1, 52) -4HCNTL, SPCNTL (1,52) -1.0, 7.0, 4HSTAP, 8.0, 

5. 0. 0. .1.0.0. .40. 00, 

KONFIG (1, 54) -4HOPTV, 0, 0,15,0, SPEC (1, 54) -0.,0.,0., 1.0,0., 

0. , 0. , 0. , 1 .0, 

KONFIG (1, 55) -4HOPTV,0, 0, 20 , 0, SPEC (1 , 55) -0 . , 0 . , 0 . , 1 . 0, 0 . , 

0. , 0. , 0. , 1 .0, 

KONFIG (1, 56) -4HOPTV, 0,0, 13 , 0, SPEC (1 , 56) -0 . , 0 . , 0 . ,1.0,0., 

0. , 0. , 0. , 1 .0, 

KONFIG (1, 58) -4HCNTL, SPCNTL (1,58) -4.0, 6. 0,4HDOUT, 6.0, 

3 . 0. 1 .0. 0. . 0. . 0. , 

KONFIG (1, 59) -4HDUCT, 4,0,5, 0, SPEC (1 , 59) -0 . , 

(END 

SD MACH-0. 30, ALTP-5000. , ETAR-1 . 00, MODE-1 , 

SD MACH-0 . 50, ALTP-10000 . , ETAR-1 . 00, MODE-1 , 

SD MACH-0 . 80, ALTP-20000 . , ETAR-1 . 00, MODE-1 , 

SD MACH-1 . 00, ALTP-30000 . , ETAR-1 . 00, MODE-1 , 

SD MACH-1 . 40, ALTP-40000 . , ETAR-0 . 9782, MODE-1 , 

SD MACH-1 . 80 , ALTP-50000 . , ETAR-0 . 9445, MODE-1 , 

SD MACH-0 . 00, ALTP-0 . , ETAR-1 . 00,MODE-2, 

XNUM-4HINTI, 4HSPLI, 4HFFNI , 4HVLVI, 4H2FNI , 4H2FNO, 4HHPCI, 4HHPCO, 4HHPTI , 
4HLPTI, 4HLPTO, 4H , 4HNOZI, 4HNOZO, 4HFFNO, 4HDNZI , 4HDNZO, 4HCOOL, 4HBDUC, 
4HBNZI, 4HBNZO, 

SPEC (4, 6) -2860, 

SPEC (7, 15) -0, SPEC (9, 51) -0, SPEC (9, 39) -0, IDONE (7) -0, 

SD MACH-0. 30, ALTP-5000. , ETAR-1. 00, MODE-2, 

SPEC (9, 51) -1 , SPEC (9, 39) -1 , SPEC (9, 33) -1, SEND 
SD MACH-0. 50, ALTP-10000., ETAR-1. 00, MODE-2, SEND 
SD MACH-0. 80, ALTP-20000. , ETAR-1. 00, MODE-2, SEND 
SD MACH-1. 00, ALTP-30000., ETAR-1. 00, MODE-2, SEND 
SD MACH-1. 40, ALTP-40000., ETAR-0. 9782, MODE-2, SEND 
SD MACH-1. 80, ALTP-50000., ETAR-0. 9445, MODE-2, SEND 


SEND 
SEND 
SEND 
SEND 
SEND 
SEND 


SEND 
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Figure 14. Block Diagram of Tandem Fan Engine, Mode 1, In Series, Mixed-Flow 
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Figure 15. Block Diagram of Tandem Fan Engine, Mode 2, Parallel, High Bypass 
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TEST C&SE: MIXED PLOW TURBOFAN SANFORD NEW 

LONG-T, NCODE— 1, DOUTHD-T, DRAW-T , 4 END 

TABLE DATA INPUT SUMMARY 10 TABLES 


12.6 EXAMP LE OUTPUT FOR SAMPT.f. tasr a 


The following is a sample output from the NNEP89 program. This output is from 
Sample 3, a simple mixed-flow, 1 spool turbofan (Section 12.3). 
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12.7 EXAMPLE COMPRES SOR MAP TABLES 


The following is a sample compressor map table with variable stator angles. Two stator 
angles (0.0, 10.0) are given in the table. The corresponding output graph of this table which 
results from setting MAPLOT=T is shown in Figure 16. The corrected speed values are labeled 
along the top of the graph and the R values are labeled on the right side of the graph. The 
efficiency contours are numbered and the corresponding values are listed in the lower right 
hand corner. 


2001 HIGH PRESSURE COMPRESSOR WITH VARIABLE STATORS 


ANGL 

2 

0.00 

10.00 






SPED 

15 

0.600 

0.700 

0.750 

0.800 

0.810 

0.820 

G.830 



0.840 

0.850 

0.860 

0.870 

0.900 

0.935 

0.985 



1.035 







R 

7 

1.000 

1.050 

1.150 

1.300 

1.450 

1.600 

1.750 

FLOW 

7 

0.3520 

0.3580 

0.3640 

0.3730 

0.3820 

0.3840 

0.3840 

FLOW 

7 

0.3910 

0.3960 

0.4C60 

0.4140 

0.4210 

0.4250 

0.4260 

FLOW 

7 

0.4330 

0. 4440 

0.4500 

0.4550 

0.4580 

0.4600 

0.4610 

FLOW 

7 

0.4690 

0.4840 

0.4930 

0.5000 

0.5040 

0.5050 

0.5060 

FLOW 

7 

0.5080 

0.5230 

0.5350 

0.5480 

0.5530 

0.5550 

0.5560 

FLOW 

7 

0.5690 

0.5820 

0.5930 

0.6080 

0.6170 

0.6210 

0.6240 

FLOW 

7 

0.6240 

0. 6380 

0.6540 

0.6700 

0.6770 

0.6800 

0.6810 

FLOW 

7 

0.6580 

0.6720 

0.6860 

0.7020 

0.7100 

0.7140 

0.7160 

FLOW 

7 

0.6880 

0.7030 

0.7160 

0.7300 

0.7370 

0.7410 

0.7440 

FLOW 

7 

0.7240 

0.7350 

0.7510 

0.7660 

0.7750 

0.7780 

0.7800 

FLOW 

7 

0.7580 

0.7710 

0.7850 

0.8020 

0.8090 

0.8110 

0.8150 

FLOW 

7 

0.8430 

0.8550 

0.8680 

0.8800 

0.8850 

0.8890 

0.8910 

FLOW 

7 

0.9120 

0.9250 

0. 9350 

0.9480 

0.9510 

0.9530 

0.9540 

FLOW 

7 

0.9860 

0.9950 

1.0040 

1.0090 

1 . 0100 

1 .0100 

1.0100 

FLOW 

7 

1.0600 

1. 0600 

1.0600 

1.0600 

1.0600 

1.0600 

1.0600 

SPED 

15 

0.600 

0.700 

0.750 

0.800 

0.810 

0.820 

0.830 



0.840 

0.850 

0.860 

0.870 

0.900 

0.935 

0.985 



1.035 







R 

7 

1.000 

1 .050 

1 .150 

1 . 300 

1 .450 

1 . 600 

1 .750 

FLOW 

n 

0.7520 

0.7580 

0.7640 

0.7730 

0.7820 

0.7840 

0.7840 

FLOW 

1 

0.7910 

0.^960 

0.8060 

0,8140 

0.8210 

0.8250 

0.8260 

FLOW 

1 

0.8330 

0.8440 

0.8500 

0.8550 

0.8380 

0.8600 

0.8610 

FLOW 

1 

0.8690 

0.8840 

0.8930 

0.9000 

0 . 9040 

0.9050 

0.9060 

FLOW 

1 

0.9080 

0.9230 

0.9350 

0.9480 

0.9530 

0.9550 

0.9560 

FLOW 

7 

0.9690 

0.9820 

0.9930 

1.0080 

1.0170 

1 .0210 

1.0240 

FLOW 

7 

1.0240 

1.0380 

1.0540 

1.0700 

1 . 0770 

1.0800 

1.0810 

FLOW 

7 

1.0580 

1.0720 

1.0860 

1.1020 

1.1100 

1.1140 

1.1160 

FLOW 

7 

1.0880 

1.1030 

1.1160 

1.1300 

1.1370 

1.1410 

1.1440 

FLOW 

7 

1.1240 

1.1350 

1.1510 

1.1660 

1.1750 

1.1780 

1.1800 

FLOW 

7 

1.1580 

1.1710 

1.1850 

1.2020 

1.2090 

1.2110 

1.2150 

FLOW 

7 

1.2430 

1.2550 

1.2680 

1.2800 

1.2850 

1.2890 

1.2910 

FLOW 

7 

1.3120 

1.2250 

1.3350 

1.3480 

1.3510 

1.3530 

1.3540 

FLOW 

7 

1.3860 

1.3950 

1.4040 

1.4090 

1.4100 

1 .4iC0 

1.4100 

FLOW 

7 

1.4600 

1.4600 

1.4600 

1.4600 

1.4600 

1.4600 

1.4600 

EOT 









2002 

HIGH PRESSURE 

COMPRESSOR WITH VARIABLE STATORS 



ANGL 

2 

0.00 

10.00 






SPED 

15 

0.600 

0.700 

0.750 

0.800 

0.810 

0.820 

0.830 



0.840 

0.850 

0.860 

0.870 

0.900 

0.935 

0.985 



1.035 







R 

7 

1.000 

1.050 

1.1L0 

1.300 

1.450 

1.600 

1.750 

EFF 

7 

0.9000 

0.8500 

0.7160 

0.5450 

0.3400 

0.2000 

0.2000 

EFF 

7 

0.9400 

0.9080 

0.8350 

0.6600 

0.4450 

0.3150 

0.2800 

EFF 

7 

0.9540 

0.9400 

0.8830 

0.7650 

0.5350 

0.2850 

0.1800 

EFF 

7 

0.9640 

0.9540 

0.9160 

0.8120 

0.6450 

0.4000 

0.2400 

EFF 

7 

0.9730 

0 9660 

0.9430 

0.8650 

0.7450 

0.5620 

0.3200 

EFF 

7 

0.9860 

0.9790 

0.9650 

0.9050 

0.8180 

0.6950 

C . 5470 

EFF 

7 

0.9960 

0.9920 

0.9820 

0.9380 

0.8680 

0.7660 

0.6350 

EFF 

7 

1.0030 

0.9980 

0.9900 

0.9500 

0.8900 

0.7890 

0.6660 

EFF 

7 

1.0070 

1.0040 

0.9970 

0.0520 

0.9030 

0.8070 

0.6810 

EFF 

7 

1.0110 

1.0090 

1.0030 

0.9710 

0.9180 

0.8300 

0.7040 
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EFF 

7 

1.0140 

1.0120 

1.0080 

0.9830 

0.9310 

0.8360 

0.7100 

EFF 

7 

1.0180 

1.0150 

1.0150 

0.9950 

0.9420 

0.8550 

0.7250 

FFF 

7 

1.0150 

1.0140 

1.0110 

0.9820 

0.9300 

0.8430 

0.7000 

EFF 

7 

1.0070 

1.0010 

0.9930 

0.9570 

0.9060 

0.8160 

0.6680 

EFF 

7 

0.9180 

G . 9180 

0.9090 

0.8900 

0.8530 

0.7740 

0.6070 

SPED 

15 

0.600 

0.840 

1.035 

0.700 

0.850 

0.750 

0.860 

0.800 

0.870 

0.810 

0.900 

0.820 

0.935 

0.830 

0.985 

R 

7 

1.000 

1.050 

1.150 

1.300 

1.450 

1.600 

1.750 

EFF 

7 

0.8550 

0.8075 

0.6802 

0.5177 

0.3230 

0.1900 

0.1900 

EFF 

7 

0.8930 

0.8626 

0.7932 

0.6270 

0.4227 

0.2993 

0.2660 

EFF 

7 

0.9063 

0.8930 

0.8388 

0.7268 

0.5082 

0.2707 

0.1710 

EFF 

7 

0.9158 

0.9063 

0 8702 

0.7714 

0.6127 

0.3800 

0.2280 

EFF 

7 

0.9244 

0.9177 

0.8959 

0.8217 

0.7077 

0.5339 

0.3040 

EFF 

7 

0.9367 

0.9301 

0 . i*l 68 

0.8597 

0.7771 

0.6603 

0.5196 

EFF 

7 

0.9462 

0.9424 

0. ^29 

0.8911 

0.8246 

0.7277 

0.6032 

EFF 

7 

0.9528 

0.9481 

0.9405 

0.9025 

0.8455 

0.7495 

0.6327 

EFF 

7 

0.9567 

0.9538 

0.9471 

0.9139 

0.8579 

0.7666 

0.6470 

EFF 

7 

0.9604 

0.9585 

0.9528 

0.9224 

0.8721 

0.7885 

0.6688 

EFF 

7 

0.9633 

0.9614 

0.9576 

0.9338 

0.8845 

0.7942 

0.6745 

EFF 

7 

0.9671 

0.9642 

0.9642 

0.9452 

0.8949 

0.8122 

0.6887 

EFF 

7 

0.9642 

0.9633 

0.9604 

0.9329 

0.8835 

0.8009 

0.6650 

EFF 

7 

0.9567 

0.9509 

0.9434 

0.9092 

0.8607 

0.7752 

0.6346 

EFF 

EOT 

1003 

ANGL 

7 

2 

0.8721 0.8721 0.8636 0.8455 0.8104 

HIGH PRESSURE COMPRESSOR WITH VARIABLE STATORS 
0.00 10.00 

0.7353 

0.5766 

SPED 

15 

0.600 

0.840 

1.035 

0.700 

0.850 

0.750 

0.860 

0.800 

0.870 

0.810 

0.900 

0.820 

0.935 

0.830 

0.985 

R 

7 

1.000 

1.050 

1.150 

1.300 

1.450 

1.600 

1 . 750 

PR 

7 

2.0730 

2.0120 

1.7360 

1 .4600 

1.1530 

1.0000 

1.0000 

PR 

7 

2.5940 

2.4410 

2.1880 

1.8280 

1.4750 

1.2840 

1.2300 

PR 

7 

3.1610 

2.9470 

2.6250 

2.2260 

1.7970 

1.5360 

1.4220 

PR 

7 

3.6060 

3.3910 

3.0650 

2.6090 

2 .1650 

1.8430 

1 . 6670 

PR 

7 

4.0810 

3.8200 

3.5290 

3.0690 

2.5940 

2.1880 

1 . 9430 

PR 

7 

4.7400 

4.4570 

4.1420 

3.6210 

3.1080 

2.6480 

2.3030 

PR 

7 

5.3230 

5.0700 

4.7400 

4.1960 

3.6060 

3.0690 

2 . 6320 

PR 

7 

5.6750 

5.4300 

5.0620 

4 . 4870 

3.8890 

3.2990 

2.8240 

PR 

7 

6.0050 

5.7670 

5.3840 

4.7550 

4.1040 

3.4910 

2 . 9620 

PR 

7 

6.4030 

6.1350 

5.7360 

5.0770 

4.4030 

3.7590 

3.1690 

PR 

7 

6.8020 

6.5030 

6.0970 

5.4070 

4.6940 

3.9660 

3.3610 

PR 

7 

7.8210 

7.4230 

6.9550 

6.1350 

5.3300 

4.5100 

3.7740 

PR 

7 

8.6490 

8.1810 

7.6600 

6.7480 

5.8510 

4 . 9700 

4.1040 

PR 

7 

9.3690 

8.9250 

8.3580 

7.3460 

6.3500 

5.3530 

4.4340 

PR 

7 

9.8140 

9.6530 

8.9630 

7.8130 

6.7710 

5.7060 

4 . 7020 

SPED 

15 

0.600 

0.840 

1.035 

0.700 

0.850 

0.750 

0.860 

0.800 

0.870 

0.810 

0.900 

0.820 

0.935 

0.830 

0.985 

R 

7 

1.000 

1.050 

1.150 

1.300 

1.450 

1.600 

1 .750 

PR 

7 

2.7919 

2.6900 

2.2291 

1.7682 

1.2555 

1.0000 

1.0000 

PR 

7 

3.6620 

3.4065 

^.9840 

2.3828 

1.7932 

1.4743 

1.3841 

PR 

7 

4.6089 

4.2515 

3.7137 

3.0474 

2.3310 

1.8951 

1.7047 

PR 

7 

5.3520 

4.9930 

4.4552 

3.6870 

2.9455 

2.4078 

2 .1139 

PR 

7 

6.1453 

5.7094 

5.2234 

4.4552 

3.6620 

2.9840 

2.5748 

PR 

7 

7.2458 

6.7732 

6.2471 

5.3771 

4.5204 

3.7522 

3.1760 

PR 

7 

8.2194 

7.7969 

7.2458 

6.3373 

5.3520 

4.4552 

3.7254 

PR 

7 

8.8072 

8.3981 

7.7835 

6.8233 

5.8246 

4.8393 

4.0461 

PR 

7 

9.3583 

8.9609 

8.3213 

7.2708 

6.1837 

5.1600 

4.2765 

PR 

7 

10.0230 

9.5754 

8.9091 

7.8086 

6.6830 

5.6075 

4.6222 

PR 

7 

10.6893 

10.1900 

9.5120 

8.3597 

7.1690 

5.9532 

4 . 9429 

PR 

7 

12.3911 

11.7264 

10.9448 

9.5754 

8.2311 

6.8617 

5.6326 

PR 

7 

13.7738 

12.9923 

12.1222 

10.5992 

9.1012 

7.6299 

6.1837 

PR 

7 

14.9762 

14.2347 

13.2879 

11.5978 

9.9345 

8.2695 

6.7348 

PR 

EOT 

7 

15.7194 

15.4505 

14.2982 

12.3777 

10.6376 

8.8590 

7.1823 
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12.8 EXAMPLE TURBINE MAP TABLES 


The following is a sample turbine map table. The corresponding output graph of this 
table which results from setting MAPLOT=T is shown in Figure 17. The corrected speed values 
are labeled with symbols and the corresponding speed values are listed in the lower right hand 
corner of the graph. 


2004 HIGH PRESSURE TURBINE 


AREA 

1 

1.00 







SPED 

3 

4000. 

5000. 

6000. 





PR 

14 

1.000 

1.300 

1.500 

1.600 

1.800 

2.000 

2.200 



2.500 

2.800 

3.100 

3.300 

3.500 

3.600 

5.000 

FLOW 

14 

0.000 

15.300 

17.100 

17.775 

18.625 

19.150 

19.460 



19.750 

19.900 

19.980 

20.010 

20.040 

20.040 

20.041 

FLOW 

14 

0.000 

15.775 

17.100 

17.575 

18.225 

18.700 

19.040 



19.360 

19.540 

19.640 

19.670 

19.700 

19.700 

19.701 

FLOW 

14 

0.000 

16.225 

17.125 

17.500 

18.040 

18.450 

18.750 

EOT 


19.050 

19.1' 0 

19.280 

19.310 

19.340 

19.340 

19.341 

2005 

HIGH PRESSURE TURBINE 






AREA 

1 

1.00 







SPED 

4 

4000.0 

5000.0 

6000.0 

8000.0 




PR 

14 

1.000 

1.250 

1.750 

2.000 

2 .150 

2.380 

2.500 



2 . 750 

3.250 

3.500 

4.000 

4.500 

4 . 750 

5 . 000 

EFF 

14 

0.8370 

0.8419 

0.8512 

0.8557 

0.8581 

0.8615 

0.8635 



0.8672 

0.8734 

0.8760 

0.8786 

0.8790 

0.8782 

0.8770 

EFF 

14 

0.8400 

0.8495 

0.8657 

0.8725 

0.8765 

0 . 8806 

0. 8811 



0.8815 

0.8819 

0.8820 

0.8802 

0.8762 

0.8733 

0 . 8700 

EFF 

14 

0.8400 

0.8492 

0.8648 

0.8705 

0.8735 

0.8772 

0.8792 



0 . 8822 

0.8867 

0.8881 

0.8891 

0.8877 

0. 8862 

0 . 8840 

EFF 

14 

0.8400 

0 . 8489 

0.8636 

0.8687 

0.8711 

0.8726 

0. 8732 

EOT 


0 .8736 

0.8727 

0.8720 

0.8705 

0.8688 

0 .8678 

0.8668 
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